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The current status of the recent developments of the second-order gauge-invariant cosmological 
perturbation theory is reviewed. To show the essence of this perturbation theory, we concentrate 
only on the universe filled with a single scalar field. Through this review, we point out the problems 
which should be clarified for the further theoretical sophistication of this perturbation theory. We 
also expect that this theoretical sophistication will be also useful to discuss the theoretical predictions 
of Non-Gaussianity in CMB and comparison with observations. 
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I. INTRODUCTION 

The general relativistic cosmological linear perturba- 
tion theory has been developed to a high degree of so- 
phistication during the last 30 years [ll-Q. One of the mo- 
tivations of this development was to clarify the relation 
between the scenarios of the early universe and cosmo- 
logical data, such as the cosmic microwave background 
(CMB) anisotropics. Recently, the first-order approxi- 
mation of our universe from a homogeneous isotropic one 
was revealed through the observation of the CMB by the 
Wilkinson Microwave Anisotropy Probe (WMAP)[3, 
the cosmological parameters are accurately measured, we 
have obtained the standard cosmological model, and the 
so-called "precision cosmology" has begun. These de- 
velopments in observations were also supported by the 
theoretical sophistication of the linear order cosmologi- 
cal perturbation theory. 

The observational results of CMB also suggest that 
the fluctuations of our universe are adiabatic and Gaus- 
sian at least in the first-order approximation. We are 
now on the stage to discuss the deviation from this first- 
order approximation from the observational [5| and theo- 
retical sides[^ 01 through the non-Gaussianity, the non- 
adiabaticity, and so on. These will be goals of future 
satellite missions. With the increase of precision of the 
CMB data, the study of relativistic cosmological pertur- 
bations beyond linear order is a topical subject. The 
second-order cosmological perturbation theory is one of 
such perturbation theories beyond linear order. 

Although the second-order perturbation theory in gen- 
eral relativity is an old topic, a general framework of 
the gauge-invariant formulation of the general relativistic 
second-order perturbation has been proposed [1, This 
general formulation is an extension of the works of Bruni 
et al.[l3l and has also been applied to cosmological per- 
turbations: The derivation of the second-order Einstein 
equation in a gauge-invariant manner without any gauge 
fixing [TT|: Applicability in more generic situations ^12J: 
Confirmation of the consistency between all components 
of the second-order Einstein equations and equations of 
motions [Tsj. We also note that the radiation case has 
recently been discussed by treating the Boltzmann equa- 
tion up to second order [1^ along the gauge-invariant 



manner of the above series of papers by the present au- 
thor. 

In this review article, we summarize the current status 
of this development of the second-order gauge-invariant 
cosmological perturbation theory through the simple sys- 
tem of a sclar field. Through this review, we point out 
the problems which should be clarified and directions of 
the further development of the theoretical sophistication 
of the general relativistic higher-order perturbation the- 
ory, especially in cosmological perturbations. We expect 
that this sophistication will be also useful to discuss the 
theoretical predictions of Non-Gaussianity in CMB and 
comparison with observations. 

The organization of this paper is as follows. In 
Sec. ini we review the general framework of the second- 
order gauge invariant perturbation theory developed in 
Refs. @, [3, [m, • This review also includes additional 
explanation not given in those papers. In Sec. IIIIl we 
also the derivations of the second-order perturbation of 
the Einstein equation and the energy-momentum tensor 
from general point of view. For simplicity, in this pa- 
per, we only consider a single scalar field as a matter 
content. The ingredients of Sec. HIl and IIIIl will be appli- 
cable to perturbation theory in any theory with general 
covariance, if the decomposition formula (|2.23p for the 
linear-order metric perturbation is correct. In Sec. lIVI we 
summarize the Einstein equations in the case of a back- 
ground homogeneous isotropic universe, which are used 
in the derivation of the first- and second-order Einstein 
equations. In Sec. El the first-order perturbation of the 
Einstein equations and the Klein-Gordon equations are 
summarized. The derivation of the second-order pertur- 
bations of the Einstein equations and the Klein- Gordon 
equations, and their consistency are reviewed in Sec. I VII 
The final section. Sec. IVIII is devoted to a summary and 
discussions. 



II. GENERAL FRAMEWORK OF THE 
GENERAL RELATIVISTIC GAUGE-INVARIANT 
PERTURBATION THEORY 

In this section, we review the general framework of 
the gau ge invariant perturbation theory developed in 
Refs. To develop the general relativistic 
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gauge-invariant perturbation theory, we first explain the 
general arguments of the Taylor expansion on a mani- 
fold without introducing an explicit coordinate system 
in Sec lII Al Further, we also have to clarify the notion 
of "gauge" in general relativity to develop the gauge- 
invariant perturbation theory from general point of view, 
which is explained in Sec. IIIBl After clarifying the no- 
tion of "gauge" in general relativistic perturbations, in 
Sec. Ill CI we explain the formulation of the general rela- 
tivistic gauge-invariant perturbation theory from general 
point of view. Although our understanding of "gauge" in 
general relativistic perturbations essentially is different 
from "degree of freedom of coordinates" as in many lit- 
erature, "a coordinate transformation" is induced by our 
understanding of "gauge" . This situation is explained in 
Sec. Ill D1 To exclude "gauge degree of freedom" which is 
unphysical degree of freedom in perturbations, we con- 
struct "gauge-invariant variables" of perturbations as re- 
viewed in Sec. Ill El These "gauge-invariant variables" 
are regarded as physical quantities. 



A. Taylor expansion of tensors on a manifold 

First, we briefly review the issues on the general form of 
the Taylor expansion of tensors on a manifold A4. The 
gauge issue of general relativistic perturbation theories 
which we will discuss is related to the coordinate transfor- 
mation. Therefore, we have to discuss the general form of 
the Taylor expansion without the explicit introduction of 
coordinate systems. Although we only consider the Tay- 
lor expansion of a scalar function / : i— )■ K, here, the 
resulting formula is extended to that for any tensor field 
on a manifold as in Appendix [X] We have to emphasize 
that the general formula of the Taylor expansion shown 
here is the starting point of our gauge-invariant formula- 
tion of the second-order general relativistic perturbation 
theory. 

The Taylor expansion of a function / is an approxi- 
mated form of /(g) at q £ in terms of the variables 
a.t p G A4, where q is in the neighborhood of p. To de- 
rive the formula for the Taylor expansion of /, we have 
to compare the values of / at the different points on 
the manifold. To accomplish this, we introduce a one- 
parameter family of diffeomorphisms ^\ : Ai A4, 
where ^x{p) — q and $a=o(p) = P- One example of a 
diffeomorphisms $a is an exponential map with a gen- 
erator. However, we consider a more general class of 
diffeomorphisms . 

The diffeomorphism induces the pull-back $^ of 
the function / and this pull-back enable us to compare 
the values of the function / at different points. Further, 
the Taylor expansion of the function f(q) is given by 
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Since this expression hold for an arbitrary smooth func- 
tion /, the function / in Eq. (|2.ip can be regarded as a 
dummy. Therefore, we may regard the Taylor expansion 
(|2.ip to be the expansion of the pull-back $^ of the diffeo- 
morphism , rather than the expansion of the function 
/• 

According to this point of view, Sonego and BrunifT^ 
showed that there exist vector fields and such that 
the expansion (|2.ip is given by 



(2.2) 



without loss of generality (see Appendix |^ . Equation 
(|2.2p is not only the representation of the Taylor ex- 
pansion of the function /, but also the definitions of 
the generators and These generators of the one- 
parameter family of diffeomorphisms $a represent the di- 
rection along which the Taylor expansion is carried out. 
The generator is the first-order approximation of the 
fiow of the diffeomorphism $a, and the generator ^2 is the 
second-order correction to this fiow. We should regard 
the generators and to be independent. Further, as 
shown in Appendix |^ the representation of the Taylor 
expansion of an arbitrary scalar function / is extended 
to that for an arbitrary tensor field Q just through the 
replacement f ^ Q- 

We must note that, in general, the representation (|2.2p 
of the Taylor expansion is different from an usual ex- 
ponential map which is generated by a vector field. In 
general. 



(2.3) 



(2.1) 



As noted in Ref. [10'|, if the second-order generator ^2 
in Eq. (|2.2I) is proportional to the first-order generator 

in Eq. (12.21) . the diffeomorphism $a is reduced to an 
exponential map. Therefore, one may reasonably doubt 
that <i>A forms a group except under very special con- 
ditions. However, we have to note that the properties 
(|2.3p does not directly mean that $a does not form a 
group. There will be possibilities that ^a form a group 
in a different sense from exponential maps, in which the 
properties (|2.3p will be maintained. 

Now, we give an intuitive explanation of the repre- 
sentation (|2.2[) of the Taylor expansion through the case 
where the scalar function / in Eq. (|2.2p is a coordinate 
function. When two points p,q G A4 in Eq. (|2.2p are 
in the neighborhood of each other, we can apply a co- 
ordinate system Ad H> R" {n — dim Ad), which denoted 
by {2;^}, to an open set which includes these two points. 
Then, we can measure the relative position of these two 
points p and g in in terms of this coordinate system 
in R" through the Taylor expansion (12. 2p . In this case, 
we may regard that the scalar function / in Eq. (|2.2I) is 
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x^'{q) ^''(q) 




FIG. 1: (a) The second term A^i(p) in Eq. (|2.4|) is the vec- 
tor which point from the point x'^{p) to the point x'^{q) in 
the sense of the first-order correction, (b) If we look at the 
neighborhood of the point x'^{q) in detail, the vector A^i(p) 
may fail to point to x'^{q) in the sense of the second order. 
Therefore, it is necessary to add the second-order correction 



a coordinate function and Eq. (|2.2p yields 

+0(A3), (2.4) 

The second term A^i (p) in the right hand side of Eq. (|2.4p 
is familiar. This is regarded as the vector which point 
from the point a;''(p) to the point x^(<?) in the sense 
of the first-order correction as shown in Figllja). How- 
ever, in the sense of the second order, this vector X£_i{p) 
may fail to point to x^{q). Therefore, it is necessary 
to add the second-order correction as shown in Fig[ljb) . 
As a correction of the second order, we may add the 
term \y? £X{p)(^v^i{p) ■ This second-order correction cor- 
responds to that comes from the exponential map which 
is generated by the vector field . However, this cor- 
rection completely determined by the vector field (J^. 
Even if we add this correction comes from the exponen- 
tial map, there is no guarantee that the corrected vector 
X^i{p)~^\X^^i{p)dv£,i{p) does point to x'^{q) in the sense 
of the second order Thus, we have to add the new cor- 
rection \}?^2{p) of ttis second order, in general. 

Of course, without this correction ^A^^2 the vector 
which comes only from the exponential map generated by 
the vector field might point to the point Ac- 
tually, this is possible if we carefully choose the vector 
field taking into account of the deviations at the sec- 
ond order. However, this means that we have to take 
care of the second-order correction when we determine 
the first-order correction. This contradicts to the phi- 
losophy of the Taylor expansion as a perturbative ex- 
pansion, in which we can determine everything order by 
order. Therefore, we should regard that the correction 
iA^^2 (p) is necessary in general situations. 



B. Gauge degree of freedom in general relativity 

Since we want to explain the gauge-invariant pertur- 
bation theory in general relativity, first of all, we have 
to explain the notion of "gauge" in general relativity [T5|. 
General relativity is a theory with general covariance, 
which intuitively states that there is no preferred coor- 
dinate system in nature. This general covariance also 
introduce the notion of "gauge" in the theory. In the 
theory with general covariance, these "gauge" give rise to 
the unphysical degree of freedom and we have to fix the 
"gauges" or to extract some invariant quantities to obtain 
physical result. Therefore, treatments of "gauges" are 
crucial in general relativity and this situation becomes 
more delicate in general relativistic perturbation theory 
as explained below. 

In 1964, Sachs 16j pointed out that there are two kinds 
of "gauges" in general relativity. Sachs called these two 
"gauges" as the first- and the second-kind of gauges, re- 
spectively. Here, we review these concepts of "gauge" . 

1. First kind gauge 

The first kind gauge is a coordinate system on a sin- 
gle manifold M.. Although this first kind gauge is not 
important in this paper, we explain this to emphasize 
the "gauge" discussing in this paper is different from this 
first kind gauge. 

In the standard text book of manifolds (for example, 
see [23} ) , the following property of a manifold is written: 
on a manifold, we can always introduce a coordinate sys- 
tem as a diffeomorphism il)a from an open set Oa C M 
to an open set i}^a{Oa) C (n = dimA^). This diffeo- 
morphism ■0c, i.e., coordinate system of the open set O^, 
is called gauge choice (of the first kind). If we consider 
another open set in C we have another gauge 
choice ■0^ : ^ tpp{Op) C M" for O^. If these two 
open sets Oa and Op have the intersection OaCiOp ^ 0, 
we can consider the diffeomorphism ippot/j^^. This diffeo- 
morphism 1^13 o is just a coordinate transformation: 

tpaiOa n O/3) C K" IppiOa H Ofj) C M" , which is 

called gauge transformation (of the first kind) in general 
relativity. 

According to the theory of a manifold, coordinate sys- 
tem are not on a manifold itself but we can always intro- 
duce a coordinate system through a map from an open set 
in the manifold to an open set of M". For this reason, 
general covariance in general relativity is automatically 
included in the premise that our spacetime is regarded 
as a single manifold. The first kind gauge does arise due 
to this general covariance. The gauge issue of the first 
kind is represented by the question, which coordinate sys- 
tem is convenient? The answer to this question depends 
on the problem which we are addressing, i.e., what we 
want to clarify. In some case, this gauge issue of the first 
kind is an important. However, in many case, it becomes 
harmless if we apply a covariant theory on the manifold. 
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2. Second kind gauge 

The second kind gauge appears in perturbation theories 
in a theory with general covariance. This notion of the 
second kind "gauge" is the main issue of this paper. To 
explain this, we have to remind what we are doing in 
perturbation theories. 

First, in any perturbation theories, we always treat two 
spacetime manifolds. One is the physical spacetime M.. 
This physical spacetime Ai is our nature itself and we 
want to describe the properties of this physical spacetime 
M through perturbations. The other is the background 
spacetime A^q. This background spacetime have noth- 
ing to do with our nature and this is a fictitious manifold 
which is prepared by us. This background spacetime is 
just a reference to carry out perturbative analyses. We 
emphasize that these two spacetime manifolds M and 
Mo are distinct. Let us denote the physical spacetime 
by [M-, gab) and the background spacetime by {M.Q,gab), 
where gab is the metric on the physical spacetime mani- 
fold, M, and gab is the metric on the background space- 
time manifold, M-a- Further, we formally denote the 
spacetime metric and the other physical tensor fields on 
by Q and its background value on A4o by Qo- 

Second, in any perturbation theories, we always write 
equations for the perturbation of the physical variable Q 
in the form 

Q(V) = Qo(p) + '5Q(p). (2.5) 

Usually, this equation is simply regarded as a relation be- 
tween the physical variable Q and its background value 
Qo, or as the definition of the deviation 5Q of the phys- 
ical variable Q from its background value Qo- However, 
Eq. (|2.5p has deeper implications. Keeping in our mind 
that we always treat two different spacetimes, AA and 
jMo, in perturbation theory, Eq. (|2.5p is a rather curi- 
ous equation in the following sense: The variable on the 
left-hand side of Eq. (|2.5p is a variable on M. , while the 
variables on the right-hand side of Eq. (12.51) are variables 
on A^o- Hence, Eq. (|2.5I) gives a relation between vari- 
ables on two different manifolds. 

Further, through Eq. (|2.5|) . we have implicitly identi- 
fied points in these two different manifolds. More specifi- 
cally, Q("p") on the left-hand side of Eq. (|2.5p is a field on 
M , and "p" G . Similarly, we should regard the back- 
ground value Qo{p) of Q{''p") and its deviation SQ{p) 
of Q ( "p" ) from Qo (p) , which are on the right-hand side 
of Eq. (|2.5p . as fields on A^o, and p S A^o- Because 
Eq. (|2.5[) is regarded as an equation for field variables, it 
implicitly states that the points "p" G M andp G Alo are 
same. This represents the implicit assumption of the ex- 
istence of a map A^o — >■ Al : p G Alo "p" G M, which 
is usually called a gauge choice (of the second kind) in 
perturbation theory |17l|. 

It is important to note that the second kind gauge 
choice between points on A^o and A4, which is estab- 
lished by such a relation as Eq. (12. 5p . is not unique to 
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FIG. 2: The second kind gauge is a point-identification be- 
tween the physical spacetime M\ and the background space- 
time Afo on the extended manifold A/". Through Eq. (|2.5p . we 
implicitly assume the existence of a point-identification map 
between Mx and Afo- However, this point-identification is 
not unique by virtue of the general covariance in the theory. 
We may chose the gauge of the second kind so that p G Mo 
and "p"G M\ is same (Xx)- We may also choose the gauge so 
that q G Alo and "p"G M\ is same (D^a). These are different 
gauge choices. The gauge transformation X\ — >■ yx is given 
by the diffeomorphism $ = X^^ o 3^;^. 

the theory with general covariance. Rather, Eq. (|2.5p in- 
volves the degree of freedom corresponding to the choice 
of the map X : A4o ^ A4. This is called the gauge de- 
gree of freedom (of the second kind) . Such a degree of 
freedom always exists in perturbations of a theory with 
general covariance. General covariance intuitively means 
that there is no preferred coordinate system in the theory 
as mentioned above. If general covariance is not imposed 
on the theory, there is a preferred coordinate system in 
the theory, and we naturally introduce this preferred co- 
ordinate system onto both A^o and A^. Then, we can 
choose the identification map X using this preferred co- 
ordinate system. However, there is no such coordinate 
system in general relativity due to the general covari- 
ance, and we have no guiding principle to choose the 
identification map X. Indeed, we may identify "p" G M 
with q G Alo {q ^ p) instead of p G Alo- In the above 
understanding of the concept of "gauge" (of the second 
kind) in general relativistic perturbation theory, a gauge 
transformation is simply a change of the map X . 

These are the basic ideas of gauge degree of freedom (of 
the second kind) in the general relativistic perturbation 
theory which are pointed out by Sacks [igland mathe- 
matically clarified by Stewart and Walker (iTf. Based on 
these ideas, higher-order perturbation theory has been 
developed in Refs. [ll, [l^-iil . 



C. Formulation of perturbation theory 

To formulate the above understanding in more detail, 
we introduce an infinitesimal parameter A for the per- 
turbation. Further, we consider the 4 -I- 1-dimensional 
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manifold Af = M x where 4 = dimM and A G R. 
The background spacetime Mq = J^\x=o ^^'^ physi- 
cal spacetime M — M\ — J^l^^x ^^^^ submanifolds 
embedded in the extended manifold Af. Each point on 
Af is identified by a pair (p, A), where p € A4\, and each 
point in Mq C A/" is identified by A = 0. 

Through this construction, the manifold TV is foliated 
by four-dimensional submanifolds A4x of each A, and 
these are difFeomorphic to A4 and Aio- The manifold TV 
has a natural differentiable structure consisting of the di- 
rect product of A4 and M. Further, the perturbed space- 
times Ai\ for each A must have the same differential 
structure with this construction. In other words, we re- 
quire that perturbations be continuous in the sense that 
A4 and A4o are connected by a continuous curve within 
the extended manifold TV. Hence, the changes of the dif- 
ferential structure resulting from the perturbation, for 
example the formation of singularities and singular per- 
turbations in the sense of fluid mechanics, are excluded 
from consideration. 

Let us consider the set of field equations 

£[Qx] = (2.6) 

on the physical spacetime A4\ for the physical variables 
Q\ on A4\. The field equation (|2.6p formally represents 
the Einstein equation for the metric on and the equa- 
tions for matter fields on Ai\. If a tensor field Q\ is 
given on each A4\, Q\ is automatically extended to a 
tensor field on TV by Q{p, A) := Qx{p), where p G A4x- 
In this extension, the field equation (|2.6p is regarded as 
an equation on the extended manifold TV. Thus, we have 
extended an arbitrary tensor field and the field equations 
(12. 6p on each A4x to those on the extended manifold TV. 

Tensor fields on TV obtained through the above con- 
struction are necessarily "tangent" to each A4\. To con- 
sider the basis of the tangent space of TV, we introduce 
the normal form and its dual, which are normal to each 
M\ in TV. These are denoted by {dX)a and (d/dX)'^, 
respectively, and they satisfy (dA)a(9/9A)° — 1. The 
form {dX)a and its dual, {d/dX)"^, are normal to any ten- 
sor field extended from the tangent space on each A4\ 
through the above construction. The set consisting of 
{dX)a, {d/dXy and the basis of the tangent space on 
each M.\\s regarded as the basis of the tangent space of 
M. 

Now, we define the perturbation of an arbitrary tensor 
field Q. We compare Q on M.\ with Qo on TWo, and it 
is necessary to identify the points oi M\ with those of 
M.Q as mentioned above. This point identification map is 
the gauge choice of the second kind as mentioned above. 
The gauge choice is made by assigning a diffeomorphism 
A'a : TV ^ TV such that Xx : Aio ^ Mx. Following the 
paper of Bruni et al. , we introduce a gauge choice Xx 
as an one- parameter groups of diffeomorphisms, i.e., an 
exponential map, for simplicity. We denote the generator 
of this exponential map by yq"' . This generator yq°' is 
decomposed by the basis on TV which are constructed 
above. Although the generator should satisfy some 



appropriate properties Q, the arbitrariness of the gauge 
choice Xx is represented by the tangential component of 
the generator to Mx- 

The pull-back X^Q^ which is induced by the expo- 
nential map A\, maps a tensor field Q on the physical 
manifold TWa to a tensor field X^Q on the background 
spacetime. In terms of this generator i)^"' , the pull-back 
X^Q is represented by the Taylor expansion 

g(r) - Q{Xx{p)) = X*xQ(p) 

- g(p) + A£^„g|p + iA2£2^Q|^ 

+0(A3), (2.7) 

where r = Xx{p) G AAx- Because p G TWq, we may 
regard the equation 

X*xQ{p) = Qoip) + A £^,Q\m„ (p) + £%Q\mo 

+0(A3) (2.8) 

as an equation on the background spacetime A4o, where 
Qo = Q\mo background value of the physical vari- 

able of Q. Once the definition of the pull-back of the 
gauge choice Xx is given, the first- and the second-order 
perturbations '^]j^Q and '^^Q of a tensor field Q under the 
gauge choice Xx are simply given by the expansion 

X*xQAm. - + A^Q + iA2(^)g + 0(a3) (2.9) 

with respect to the infinitesimal parameter A. Comparing 
Eqs. ()2.8|) and ()2.9|) . we define the first- and the second- 
order perturbations of a physical variable Qx under the 
gauge choice Xx by 

WQ:=£^,Q|^„, ^^^Q:^ £%Q\j^^. (2.10) 

We note that all variables in Eq. (j2.9|) are defined on TMq. 

Now, we consider two different gauge choices based on 
the above understanding of the second kind gauge choice. 
Suppose that Xx and yx are two exponential maps with 
the generators yq'^ and yq"^ on Af , respectively. In other 
words, Xx and yx are two gauge choices (see Fig. [2]). 
Then, the integral curves of each yq°- and yq'^ in TV are 
the orbits of the actions of the gauge choices Xx and 3^^ j 
respectively. Since we choose the generators yq"' and yq"^ 
so that these are transverse to each AAx everywhere on 
TV, the integral curves of these vector fields intersect with 
each AAx- Therefore, points lying on the same integral 
curve of either of the two are to be regarded as the same 
point within the respective gauges. When these curves 
are not identical, i.e., the tangential components to each 
A4x oi p^°- and yq"" are different, these point identifica- 
tion maps Xx and yx are regarded as two different gauge 
choices. 

We next introduce the concept of gauge invariance. In 
particular, in this paper, we consider the concept of order 
by order gauge invariance]!^. Suppose that Xx and yx 



6 



are two different gauge choices which are generated by 
the vector fields and j^", respectively. These gauge 
choices also pull back a generic tensor field Q on A/" to 
two other tensor fields, X^Q and y^Q, for any given 
value of A. In particular, on Aio, wc now have three 
tensor fields associated with a tensor field Q; one is the 
background value Qq of Q, and the other two are the 
pulled-back variables of Q from M\ to Mo by the two 
different gauge choices. 



Qo + A^i^g + -A2(^)Q + 0(A3) (2.11) 



= Qo + X%^Q + lx^%^Q + 0{X^) (2.12) 

Here, we have used Eq. (|2.9I) . Because X\ and yx are 
gauge choices which map from A^o to M.x, xQ\ and 
yQ\ are the different representations on in the two 
different gauges of the same perturbed tensor field Q on 
A^A- The quantities '^^^Q and '^yQ in Eqs. (|2.1ip and 
()2.12p are the perturbations of 0(A:) in the gauges X\ and 
y\ , respectively. We say that the A:th-order perturbation 

xQ of Q is order by order gauge invariant iff for any two 
gauges Xx and the following holds: 



(2.13) 



Now, we consider the gauge transformation rules be- 
tween different gauge choices. In general, the represen- 
tation "^Qx on Mq of the perturbed variable Q on Mx 
depends on the gauge choice Xx ■ If we employ a differ- 
ent gauge choice, the representation of Qx on A4o may 
change. Suppose that Xx and yx are different gauge 
choices, which are the point identification maps from A4q 
to Mx, and the generators of these gauge choices are 
given by and yjy", respectively. Then, the change 
of the gauge choice from Xx to yx is represented by the 
diffcomorphism 



$A --^{Xxr^oyx. 



(2.14) 



This diffcomorphism "I>a is the map $a : Mq — > Mo for 
each value of A e M. The diffcomorphism $a does change 
the point identification, as expected from the understand- 
ing of the gauge choice discussed above. Therefore, the 
diffcomorphism $a is regarded as the gauge transforma- 
tion <&x ■■ Xx ^ yx- 

The gauge transformation $a induces a pull-back from 
the representation p(Qx of the perturbed tensor field Q 
in the gauge choice Xx to the representation yQx in the 
gauge choice yx- Actually, the tensor fields p^Qx and 
yQx, which are defined on Mo, are connected by the 
linear map $^ as 



Mo 



= {X^^yxTiX^Q) ^<^IaQx- (2.15) 



According to generic arguments concerning the Taylor 
expansion of the pull-back of a tensor field on the same 
manifold, given in §11 A[ it should be possible to express 
the gauge transformation ^*^xQx in the form 



+0{\^), (2.16) 



where the vector fields and ^2 are the generators of 
the gauge transformation $a (sec Eq. (|2.2p ). 

Comparing the representation (j2.16l) of the Taylor ex- 
pansion in terms of the generators and of the pull- 
back ^*xxQ and that in terms of the generators ^fl" and 
yri" of the pull-back J^* o (A""^) * xQ (= $aa^), we read- 
ily obtain explicit expressions for the generators and 
^2 of the gauge transformation (f> = X^^ o yx in terms 
of the generators and yrj"" of each gauge choices as 
follows: 



C2 = [y^,A^]^ 



(2.l7) 



Further, because the gauge transformation <I>a is a map 
within the background spacetime Mo, the generator 
should consist of vector fields on A^o- This can be sat- 
isfied by imposing some appropriate conditions on the 
generators yif and ;e'7°- 

We can now derive the relation between the pertur- 
bations in the two different gauges. Up to second or- 
der, these relations are derived by substituting (|2.1ip and 
((2T2)) into (|2T6)) : 



(2)g _ or (1), 



(2.18) 
h- (2.19) 



Mo 



Here, we should comment on the gauge choice in the 
above explanation. We have introduced an exponential 
map Xx (or 3^a) as the gauge choice, for simplicity. How- 
ever, this simplified introduction of Xx as an exponential 
map is not essential to the gauge transformation rules 
(|2.18p and (|2.19p . Actually, we can generalize the dif- 
fcomorphism Xx from an exponential map. For example, 
the diffcomorphism whose pull-back is represented by the 
Taylor expansion (|2.2p is a candidate of the generaliza- 
tion. If we generalize the diffcomorphism Xx, the repre- 
sentation (|2.8p of the puUcd-back variable X^Q{p), the 
representations of the perturbations (|2.10p . and the rela- 
tions (|2.17p between generators of $a, Xx, and yx will be 
changed. However, the gauge transformation rules (|2.18p 
and (|2.19p are direct consequences of the generic Taylor 
expansion (|2.16p of 'I'a. Generality of the representa- 
tion of the Taylor expansion (|2.16p of ^a implies that 
the gauge transformation rules (|2.18p and (|2.19p will not 
be changed, even if we generalize the each gauge choice 
Xx- Further, the relations (|2.17p between generators also 
imply that, even if we employ simple exponential maps 
as gauge choices, both of the generators £,f and ^2 are 
naturally induced by the generators of the original gauge 
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choices. Hence, we conclude that the gauge transforma- 
tion rules (|2.18p and (|2.19p are quite general and irre- 
ducible. In this paper, we review the development of a 
second-order gauge-invariant cosmological perturbation 
theory based on the above understanding of the gauge 
degree of freedom only through the gauge transforma- 
tion rules (|2.18l) and (|2.19p . Hence, the developments 
of the cosmological perturbation theory presented below 
will not be changed even if we generalize the gauge choice 
Xx from a simple exponential map. 

We also have to emphasize the physical implication of 
the gauge transformation rules (|2.18l) and (j2.19p . Ac- 
cording to the above construction of the perturbation 
theory, gauge degree of freedom, which induces the trans- 
formation rules (I2.18P and (|2.19p . is unphysical degree of 
freedom. As emphasized above, the physical spacetime 
A4\ is our nature itself, while there is no background 
spacetime A4o in our nature. The background spacetime 
Mo is a fictitious spacetime and it have nothing to do 
with our nature. Since the gauge choice Xx just gives 
a relation between Mx and Mo, the gauge choice Xx 
also have nothing to do with our nature. On the other 
hand, any observations and experiments are carried out 
only on the physical spacetime Mx through the physical 
processes on the physical spacetime Mx- Therefore, any 
direct observables in any observations and experiments 
should be independent of the gauge choice Xx, i.e., should 
be gauge invariant. Keeping this fact in our mind, the 
gauge transformation rules (|2.18p and (|2.19p imply that 

the perturbations '"^Q and ^'^Q include unphysical de- 
gree of freedom, i.e., gauge degree of freedom, if these 
perturbations are transformed as (|2.18p or (|2.19p under 
the gauge transformation Xx — ^ 3^a ■ If the perturbations 
^l^Q and ^^Q are independent of the gauge choice, these 
variables are order by order gauge invariant. Therefore, 
order by order gauge-invariant variables does not include 
unphysical degree of freedom and should be related to 
the physics on the physical spacetime Mx- 



general relativistic gauge-invariant perturbation theory. 
However, we comment on the relations between the co- 
ordinate transformation, briefly. Details can be seen in 
Refs. [i,[ii,[23. 

To see that the gauge transformation of the second 
kind induces the coordinate transformation, we intro- 
duce the coordinate system {Oa, V'a} on the "background 
spacetime" Mo, where Oa are open sets on the back- 
ground spacetime and V'a are diffeomorphisms from Oa 
to (4 = dim A^o)- The coordinate system {Oa,ipa} is 
the set of the collection of the pair of open sets Oa and 
diffeomorphism Oa i— )• K'^. If we employ a gauge choice 
Xx, we have the correspondence of A^a and Mo- To- 
gether with the coordinate system tpa on Mo, this corre- 
spondence between Mx and Mo induces the coordinate 
system on A^a- Actually, Xx{Oa) for each a is an open 
set of A^A- Then, ^pa ° X^^ becomes a diffeomorphism 
from an open set XxiOa) C Mx to E^. This diffeomor- 
phism ipa o X^^ induces a coordinate system of an open 
set on A^A- 

When we have two different gauge choices Xx and yx, 
Ipa ° X^^ and Tpa ° yx^ become different coordinate sys- 
tems on A^A- We can also consider the coordinate trans- 
formation from the coordinate system o X^^ to an- 
other coordinate system tpa°yx^- Since the gauge trans- 
formation Xx yx is induced by the diffeomorphism $a 
defined by Eq. (|2.14p . the induced coordinate transforma- 
tion is given by 

y^{q):^x^{p)^[{^-'yx^){q) (2.20) 

in the passive point of view[l,[l^[23]- If we represent this 
coordinate transformation in terms of the Taylor expan- 
sion in Sec. Ill A[ up to third order, we have the coordinate 
transformation 

+0{X^)- (2.21) 



D. Coordinate transformations induced by the 
second kind gauge transformation 

In many literature, gauge degree of freedom is regarded 
as the degree of freedom of the coordinate transforma- 
tion. In the linear-order perturbation theory, these two 
degree of freedom are equivalent with each other. How- 
ever, in the higher order perturbations, we should regard 
that these two degree of freedom are different. Although 
the essential understanding of the gauge degree of free- 
dom (of the second kind) is as that explained above, the 
gauge transformation (of the second kind) also induces 
the infinitesimal coordinate transformation on the phys- 
ical spacetime A^a as a result. In many case, the un- 
derstanding of "gauges" in perturbations based on co- 
ordinate transformations leads mistakes. Therefore, we 
did not use any ingredient of this subsection in our se- 
ries of papers [g. [ol llTI - fTsj concerning about higher-order 



E. Gauge-invariant variables 

Here, inspecting the gauge transformation rules (|2.18p 
and (|2.19p . we define the gauge invariant variables for a 
metric perturbation and for arbitrary matter fields (ten- 
sor fields). Employing the idea of order by order gauge 
invariance for perturbations |T^ , we proposed a proce- 
dure to construct gauge invariant variables of higher- 
order perturbations [8]. This proposal is as follows. First, 
we decompose a linear-order metric perturbation into its 
gauge invariant and variant parts. The procedure for de- 
composing linear-order metric perturbations is extended 
to second-order metric perturbations, and we can decom- 
pose the second-order metric perturbation into gauge in- 
variant and variant parts. Then, we can define the gauge 
invariant variables for the first- and second-order pertur- 
bations of an arbitrary field other than the metric by us- 
ing the gauge variant parts of the first- and second-order 
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metric perturbations. Although the procedure for finding 
gauge invariant variables for linear-order metric pertur- 
bations is highly non-trivial, once we know this proce- 
dure, we can easily define the gauge invariant variables 
of a higher-order perturbation through a simple exten- 
sion of the procedure for the linear-order perturbations. 

Now, we review the above strategy to construct gauge- 
invariant variables. To consider a metric perturbation, 
we expand the metric on the physical spacetime Mx, 
which is pulled back to the background spacetime A4o 
using a gauge choice in the form given in (12.91) : 

Xxfjab = 9ab + \xhab+Y^'^b+0^{X), (2.22) 

where gab is the metric on Mq. Of course, the expan- 
sion (I2.22P of the metric depends entirely on the gauge 
choice X\. Nevertheless, henceforth, we do not explicitly 
express the index of the gauge choice Xx in an expression 
if there is no possibility of confusion. 

Our starting point to construct gauge invariant vari- 
ables is the assumption that we already know the proce- 
dure for finding gauge invariant variables for the linear 
metric perturbations. Then, a linear metric perturbation 
hab is decomposed as 

hab ='■ T-Lab + £xgab, (2.23) 

where T-lab and X°' are the gauge invariant and variant 
parts of the linear-order metric perturbations, i.e., under 
the gauge transformation (|2.18p . these are transformed 
as 

yUab - x-Hab = 0, yX'^ ~ xX'' = ^1 • (2.24) 

The first-order metric perturbation (|2.23p together with 
the gauge transformation rules (|2.24p does satisfy the 
gauge transformation rule (|2.18p for the first-order metric 
perturbation, i.e., 

''yhab ~ ''x^ab ~ ^'iidab- (2.25) 

As emphasized in our series of papers [1, [Tll - [T3{ . 
the above assumption is quite non-trivial and it is not 
simple to carry out the systematic decomposition (|2.23p 
on an arbitrary background spacetime, since this proce- 
dure depends completely on the background spacetime 
(A4o,5ati)- However, as we will show below, this pro- 
cedure exists at least in the case of cosmological per- 
turbations of a homogeneous and isotropic universe in 
Sec. TO 

Once we accept this assumption for linear-order metric 
perturbations, we can always find gauge invariant vari- 
ables for higher-order perturbations [8|. According to the 
gauge transformation rule (j2.19l) , the second-order metric 
perturbation Igh is transformed as 

^"y^ab — ^"^ab — '^£(,1 xhab + {"^^2 + -^li } 9ab (2.26) 

under the gauge transformation $a = {Xx)~^°yx ■ ^x 
yx- Although this gauge transformation rule is slightly 



complicated, inspecting this gauge transformation rule, 
we first introduce the variable Lab defined by 

Lab ■— lab — '2£xhab + £x9ab- (2.27) 

Under the gauge transformation — (-^a)^^ ° yx ■ 
Xx yx, the variable Lab is transformed as 

yLab — xLab = £agab, (2.28) 

a" := Q + iCuX]". (2.29) 

The gauge transformation rule (|2.28p is identical to that 
for a linear metric perturbation. Therefore, we may apply 
the above procedure to decompose hab into Hab and Xa 
when we decompose of the components of the variable 
Lab- Then, Lab can be decomposed as 

Lab ~ ^ab + £Ygab, (2.30) 

where Cab is the gauge invariant part of the variable Lab, 
or equivalently, of the second-order metric perturbation 
lab, and y is the gauge variant part of Lab, i-e., the gauge 
variant part of lab- Under the gauge transformation $a = 
{Xx)^^ o yx, the variables Lab and Y°- are transformed 
as 

yCab — xC-ab = 0, yYa — yYa = a a, (2-31) 

respectively. Thus, once we accept the assumption (j2.23p . 
the second-order metric perturbations are decomposed as 

lab —■ Lab + '2£ xhab + {£y ~ £x) 9ab, (2.32) 

where Cab and are the gauge invariant and variant 
parts of the second order metric perturbations, i.e., 

yCab - xCab = 0, yY- - xY'' = + Ki, ^l"- (2.33) 

Furthermore, as shown in Ref. [8|], using the first- and 
second-order gauge variant parts, X° and , of the met- 
ric perturbations, the gauge invariant variables for an ar- 
bitrary field Q other than the metric are given by 

(i)Q-£xQo, (2.34) 
(2)Q (2)g_2£x(^^Q-{i;y-£|}Qo.(2.35) 

It is straightforward to confirm that the variables 
defined by (|2.34p and (|2.35p are gauge invariant under 
the gauge transformation rules (|2.18p and (|2.19p . respec- 
tively. 

Equations (|2.34p and (|2.35p have very important im- 
plications. To see this, we represent these equations as 

Wq = «Q+£xQo, (2.36) 

(')Q = (2)Q + 2£x^'^Q+{i^y-£|}Qo. (2.37) 

These equations imply that any perturbation of first- 
and second-order can always be decomposed into gauge- 
invariant and gauge-variant parts as Eqs. (|2.36p and 
(|2.37p . respectively. These decomposition formulae 
(|2.36p and (|2.37l) are important ingredients in the gen- 
eral framework of the second-order general relativistic 
gauge-invariant perturbation theory. 
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III. PERTURBATIONS OF THE FIELD 
EQUATIONS 

In terms of the gauge invariant variables defined last 
section, we derive the field equations, i.e., Einstein equa- 
tions and the equation for a matter field. To derive the 
perturbation of the Einstein equations and the equation 
for a matter field (Klein-Gordon equation), first of all, we 
have to derive the perturbative expressions of the Ein- 
stein tensorlol]. This is reviewed in Sec. IIII Al We also 
derive the first- and the second-order perturbations of the 
energy momentum tensor for a scalar field and the Klein- 
Gordon equation [T^ in Sec. IIII Bl Finally, we consider 
the first- and the second-order the Einstein equations in 

Sec. Unci 

A. Perturbations of the Einstein curvature 

The relation between the curvatures associated with 
the metrics on the physical spacetime A4\ and the back- 
ground spacetime A4q is given by the relation between 
the pulled-back operator X^Va {^x^) of the covariant 
derivative Vq associated with the metric gat on and 
the covariant derivative Vq associated with the met- 
ric gab on Mq. The pulled-back covariant derivative 
X^Va (-^a"^) depends on the gauge choice Xx- The 
property of the derivative operator <-f^Va {X^^") as the 
covariant derivative on A^a is given by 

X^^a{{Xx')*X^gab)^0, (3.1) 

where Xxgab is the pull-back of the metric on A^a, which 
is expanded as Eq. (|2.22p . In spite of the gauge depen- 
dence of the operator XxVa {'^x^) i we simply denote 
this operator by Vq, because our calculations are carried 
out only on Mo in the same gauge choice Xx- Further, 
we denote the pulled-back metric Xxgab on A^a by (jab, 
as mentioned above. 

Since the derivative operator Vq (= X*\7a (X^^) ) 
may be regarded as a derivative operator on AJq that 
satisfies the property (|3.ip . there exists a tensor field 6"^^^^ 



on M-Q such that 

V.CJb = Va^b - C^ab^c, (3.2) 

where cj^ is an arbitrary one- form on Mq. From the 
property (j3.ip of the covariant derivative operator Vq on 
Mx, the tensor field C^at '^^ -^o is given by 

= i^a9db + "^bBda - ^d9ab) , (3.3) 

where g"^^ is the inverse of f^af, (see Appendix IB|). We note 
that the gauge dependence of the covariant derivative 
Va appears only through C^ab- The Riemann curvature 
R^ hf. on AIa, which is also pulled back to Ato, is given 
byll: 

R-abc = Rabc ~ 2V[^C"^f,]^ -I- 2C"'^[^C"^j]g, (3.4) 

where Rabc i^ the Riemann curvature on AJq- The per- 
turbative expression for the curvatures are obtained from 
the expansion of Eq. p.4p through the expansion of C^^^. 

The first- and the-second order perturbations of the 
Riemann, the Ricci, the scalar, the Weyl curvatures, and 
the Einstein tensors on the general background space- 
time are summarized in Ref. Q. We also derived the 
perturbative form of the divergence of an arbitrary ten- 
sor field of second rank to check the perturbative Bianchi 
identities in Ref. @. In this paper, we only present the 
perturbative expression for the Einstein tensor, and its 
derivations in Appendix IB] 

We expand the Einstein tensor Ga^ '■— Ra^ ~ ^^J'R 
A^A as 

Ga' = Ga' + A(i)G/ + h'^'^Ga' + 0{X'). (3.5) 

As shown in Appendix [B] each order perturbation of the 
Einstein tensor is given by 

WG,^ = ^'^ga'[n] + £xGa\ (3.6) 
+ 2£x'^'^Ga'+{£Y^£ji}Ga\ (3.7) 

where 



('^Ga' [A] := [A] - ha'^'K^ [A] , (DS/ [A] := -2V[„F,j'"' [A] A^^Rac, (3.8) 

^'^Ga'iA^B] := (2)s^;>[A,i?]-i^/(2)E^^[A,B], (3.9) 

^'^^a'iAB] := 2RadBj'A''^^ + 2H^/^A]H,^\[B]+2H^/^B]H,^\[A] 

+2A/\/[aH,^'^ [B] + 2B/V[aH,^'"^ [A] + 2A,'V[aH,f' [B] + 2B,'V[aH,f' [A] , (3.10) 

Hab'lA] ■■= V(,A,)- - ivM,b, (3.11) 

Habc [A] := g^dHab " [A] , i//^ [A] g^'^ff,/ [A] , [A] g^aHa'" [A] . (3.12) 
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We note that '■^•'Sa'' [*] and ^^■'^a^I*'*] i^- Eqs. p.6p and p.7p are the gauge invariant parts of the perturbative 
Einstein tensors, and Eqs. (I3.6P and p.7p have the same forms as Eqs. (|2.34p and (|2.37p . respectively. The expression 
of ^'^^Ga'' [^i B] in Eq. p.9p with Eq. p.lOp is derived by the consideration of the general relativistic gauge-invariant 
perturbation theory with two infinitesimal parameters in Refs. 's', 'o']. 

We also note that ^^'^Ga" [*] and (2)5^6 ^] defined by Eqs. ((5:H)) - (|3J0| satisfy the identities 

^J'^G.^iA] = -H,,-[A]G,^ + H,,^A]G,^, (3.13) 
V„ " [A, B] = -H,,^ [A] ^ [B] H,, ^ [B] = [A] + H,, ^ [A] ^'^G, ^ [B] + H,, ^ [B] ^ [A] 

- {Hbad [B] A"- + Hbad [A] B"'^) + {Head [B] A'^" + Head {A\ B'^") G,^ (3.14) 



for arbitrary tensor fields Aat and Bat, respectively. We 
can directly confirm these identities without specifying 
arbitrary tensors Aab and Bab of the second rank, re- 
spectively. This implies that our general framework 
of the second-order gauge invariant perturbation theory 
discussed here gives a self-consistent formulation of the 
second-order perturbation theory. These identities p.l3p 
and p.l4|) guarantee the first- and second-order pertur- 
bations of the Bianchi identity V^Gj' = and are also 
useful when we check whether the derived components of 
Eqs. dUl) and dSJ]) are correct. 

B. Perturbations of the energy momentum tensor 
and Klein-Gordon equation 

Here, we consider the perturbations of the energy mo- 
mentum tensor of the equation of motion. As a model 
of the matter field, we only consider the scalar field, for 
simplicity. Then, equation of motion for a scalar field is 
the Klein- Gordon equation. 

The energy momentum tensor for a scalar field (p is 
given by 

T/ = V„^VV - Is^' (Ve^V^^ + 2V{^)) , (3.15) 

where V{(p) is the potential of the scalar field f>. We 
expand the scalar field if as 

<f^ip + Xipi + ^X'^ip2 + 0{X^), (3.16) 



where (p is the background value of the scalar field if. 
Further, following to the decomposition formulae (j2.34p 
and (|2.35p . each order perturbation of the scalar field if 
is decomposed as 

ipi =: ipi + £xV'i (3-17) 
02 ^2 + 2£x0i + {£y - £x)^, (3.18) 

where tpi and ip2 are the first- and the second-order 
gauge-invariant perturbations of the scalar field, respec- 
tively. 

Through the perturbative expansions p.l6p and (jB2p 
of the scalar field ip and the inverse metric, the energy 
momentum tensor (j3.15p is also expanded as 

T/ = + A(i) (r, ") + i A2(2) (T„ + 0(A3). (3.19) 

The background energy momentum tensor is given by 
the replacement — > in Eq. p.isp . Further, through 
the decompositions ((O^ . ((XTTl) . and (PT5|) . the 

perturbations of the energy momentum tensor '-^^ {T^ 
and '•^^Tj) are also decomposed as 

(i)(T/) (i)r/ + i;xT/, (3.20) 

+ (£y-£|)T/, (3.21) 

where the gauge-invariant parts ^ and of the 

first and the second order are given by 



(ih- b 

' a 
(2h- b 

' n 



(3.22) 



bd^ 



1. 



-|-Vc(/5l V((5l -I- (p2 



dV 

dip 



(3.23) 



We note that Eq. (|3.20l) and (|3.21l) have the same form as (|2.36p and (|2.37l) . respectively. 
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Next, we consider the perturbation of the Klein- 
Gordon equation 

:=V"Va(^-|^(^)=0. (3.24) 

Through the perturbative expansions (|3.16p and (|2.22p . 
the Klein-Gordon equation p.24p is expanded as 

CiK) =■■ C^K) + A c'l) +^\^ +0(A3). (3.25) 



C(x) is the background Klein-Gordon equation 



(3.26) 



(1) 

The first- and the second-order perturbations C(^x) ^nd 

(2) 

C(x) also decomposed into the gauge-invariant and 
the gauge- variant parts as 



J 



where 



(1) (1) (2) (2) / 2\ 

C(K)^-C(K) +£xC(^K), C(K)^-C(K) +'^£x C(^K) + [£y - £x} C{K), (3.27) 



c'l) ■■= V"V,(pi-i?/^[?^]V,^--«"''V,Vfc(/.-^i^(^), (3.28) 

-/:°VaVfc(^ + 2-H\n''^Va^b'P - 27^°VaV6(^l - - (3.29) 



Here, we note that Eqs. p.27p have the same form as 
Eqs. ([236]) and ((237| . 

By virtue of the order by order evaluations of the Klein- 
Gordon equation, the first- and the second-order pertur- 
bation of the Klein-Gordon equation are necessarily given 
in gauge-invariant form as 

(1) (2) 

C(K)^ 0, C(K)= 0. (3.30) 

We should note that, in Ref. [l^j, we summarized the 
formulae of the energy momentum tensors for an perfect 
fluid, an imperfect fluid, and a scalar field. Further, we 
also summarized the equations of motion of these three 
matter fields: i.e., the energy continuity equation and 
the Euler equation for a perfect fluid; the energy con- 
tinuity equation and the Navier-Stokcs equation for an 
imperfect fluid; the Klein-Gordon equation for a scalar 
field. All these formulae also have the same form as the 
decomposition formulae (j2.36p and (|2.37p . In this sense, 
we may say that the decomposition formulae p.36p and 
(|2.37p are universal. 



C. Perturbations of the Einstein equation 

Finally, we impose the perturbed Einstein equation of 
each order, 

(i)G^'' = 8^G (i)r/, ^^'^Gj' = ^TiG (2)r/. (3.31) 



Then, the perturbative Einstein equation is given by 

^^^Ga[H] = SttG^i)?;'' (3.32) 

at linear order and 

(^^^a^i^l + ^'^eaM^,^] = 8^G (2)7-/ (3.33) 

at second order. These explicitly show that, order by or- 
der, the Einstein equations are necessarily given in terms 
of gauge invariant variables only. 

Together with Eqs. (I3.30p . we have seen that the 
first- and the second-order perturbations of the Einstein 
equations and the Klein-Gordon equation are necessarily 
given in gauge-invariant form. This implies that we do 
not have to consider the gauge degree of freedom, at least 
in the level where we concentrate only on the equations 
of the system. 

We have reviewed the general outline of the second- 
order gauge invariant perturbation theory. We also note 
that the ingredients of this section are independent of 
the explicit form of the background metric gab , except for 
the decomposition assumption (|2.23p for the linear-order 
metric perturbations and are valid not only in cosmologi- 
cal perturbation case but also the other generic situations 
if Eq. (I2.23|) is correct. Within this general framework, 
we develop a second-order cosmological perturbation the- 
ory in terms of the gauge invariant variables. 
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IV. COSMOLOGICAL BACKGROUND 
SPACETIME AND EQUATIONS 



EQUATIONS FOR THE FIRST-ORDER 
COSMOLOGICAL PERTURBATIONS 



The background spacetime Mq considered in cosmo- 
logical perturbation theory is a homogeneous, isotropic 
universe that is fohated by the three-dimensional hyper- 
surface S(?7), which is parametrized by rj. Each hyper- 
surface of S(?7) is a maximaUy symmetric three-space[25j, 
and the spacetime metric of this universe is given by 

gab = a^iv) {-{dv)aidr])b + jij{dx'')aidx^)b) , (4.1) 

where a = a{r]) is the scale factor, jij is the metric on the 
maximally symmetric 3-space with curvature constant 
and the indices i, j, k, ... for the spatial components run 
from 1 to 3. 

To study the Einstein equation for this background 
spacetime, we introduce the energy-momentum tensor for 
a scalar field, which is given by 



T 



Va<pVV-^'5„.'(Vc(^VV + 2yM) (4.2) 



drj 



(4.3) 



where we assumed that the scalar field ip is homogeneous 
= <fi{v) (4.4) 

and 7^'' are defined as 

lab l^J{dx')a{dx^)b, la' 1, ' {dx')a{d / Ox^ )" .{4.5) 

The background Einstein equations Gj' — ^ttGTJ' for 
this background spacetime filled with the single scalar 
field are given by 



W +K = 



V{v) j , (4.6) 
2d,^n + n^+K = 87rG(^-^{dr,(p)^ + a^V{ip)^ . (4.7) 



We also note that the equations (|4.6I) and (14.71) lead to 



(4.8) 

Equation (|4.8p is also useful when we derive the pertur- 
bative Einstein equations. 

Next, we consider the background Klein- Gordon equa- 
tion which is the equation of motion Va?], " = for the 
scalar field 



dV 

dlip + 2'Hd,M + a^— 
' dip 



(4.9) 



The Klein- Gordon equation (|4.9p is also derived from the 
Einstein equations (14. 6p and (|4.7p . This is a well known 
fact and is just due to the Bianchi identity of the back- 
ground spacetime. However, these types of relation are 
useful to check whether the derived system of equations 
is consistent. 



On the cosmological background spacetime in the last 
section, we develop the perturbation theory in the gauge- 
invariant manner. In this section, we summarize the 
first-order perturbation of the Einstein equation and the 
Klein-Gordon equations. In Sec. IV A[ we show that the 
assumption on the decomposition (j2.23p of the linear- 
order metric perturbation is correct. In Sec. IV Bl we 
summarize the first-order perturbation of the Einstein 
equation. Finally, in Sec. IV CI we show the first-order 
perturbation of the Klein-Gordon equation. 



A. Gauge-invariant metric perturbations 

Here, we consider the first-order metric perturbation 
hah and show the assumption on the decomposition (j2.23p 
is correct in the background metric Eq. (|4.1[) . To ac- 
complish the decomposition (I2.23p . first, we assume the 
existence of the Green functions A^^ {D^Di)~^, 
(A + 2K)^^, and (A + 3K)^^, where Di is the covariant 
derivative associated with the metric lij and K is the cur- 
vature constant of the maximally symmetric three space. 
Next, we consider the decomposition of the linear-order 
metric perturbation hab as 

hab = hrjrjidri)a{dr])b 

+2 (A/i(vL) + h(^v)^) {dri)ia{dx% (5.1) 



+ 2D(,h(TV)]) + h(TT)ij] {dx')a{dx^)b, 
where h(y)i, h(^TV)j7 and h(TT)ij satisfy the properties 

D'h^v), = 0, D'h(TV)^ = 0, 

h{TT)ij = hTT)j^, hT) , — ^'^/^(T)^^- = 0, (5.2) 

D''h(TT)ij = 0. 

The gauge-transformation rules for the variables /i^^, 
h{vL), fe(y )t, fe( L), h(TL), h(TV)j and h(TT)ij are derived 
from Eq. (|2.25p . Inspecting these gauge-transformation 
rules, we define the gauge- variant part Xa in Eq. (|2.23l) : 

Xa ■■= {h(VL) - ^a^d.rjh(^TL)^ {dv)a 

W (^h^TV)^ + ^A/l(TL)) {dx')a. (5.3) 

We can easily check this vector field Xa satisfies 
Eq. (12.241) . Subtracting gauge variant-part £x9ab from 
hab^ we have the gauge- invariant part T-Lab in Eq. (|2.23p : 



(1) 



(1) 



Uab = a^\-2 ^ {dr^)a{dri)b + 2 v, {dr^)(a{dx'')i,) 



(1) (1) 
-2 7y+ X^j ) {dx')a{dx^)b \ , (5.4) 
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where the properties i ■— 7*-'-Di j — Xi^ '■— 7'"' Xij 
. (1) 

= X ij= are satisfied as consequences of Eqs. (|5.2I) . 

Thus, we may say that our assumption for the de- 
composition p.23p in linear-order metric perturbation is 
correct in the case of cosmological perturbations. How- 
ever, we have to note that to accomphsh Eq. (|2.23p . 
we assumed the existence of the Green functions A~^, 
(A + 2K)-^, and (A + 3K)~\ As shown in Ref. [ul, 
this assumption is necessary to guarantee the one to 
one correspondence between the variables {hjj^,hi^,hij} 

and {/l,,,,, /l(yL), /l(TL); ^(TV)ji 'l(TT)ij}j but 

excludes some perturbative modes of the metric pertur- 
bations which belong to the kernel of the operator A, 
(A -I- 2K), and (A + 3K) from our consideration. For 
example, homogeneous modes belong to the kernel of the 
operator A and are excluded from our consideration. If 
we have to treat these modes, the separate treatments 
are necessary. In this paper, we ignore these modes, for 
simplicity. 

We also note the fact that the definition (I2.23P of 
the gauge-invariant variables is not unique. This comes 
from the fact that we can always construct new gauge- 
invariant quantities by the combination of the gauge- 
invariant variables. For example, using the gauge- 
invariant variables <I> and j/^ of the first-order met- 
ric perturbation, we can define a vector field Za by 

Za ~a $ {dr])a + a W {dx'')a which is gauge-invariant. 
Then, we can rewrite the decomposition formula (|2.23p 
for the linear-order metric perturbation as 

hab = 'Hab — £zgab + £zgab + £xgab, 

= : ICab + £x+zgab, (5.5) 

where we have defined new gauge-invariant variable ICab 
by ICab := Hab ~ £zgab- Clearly, ICab is gauge - invariant 
and the vector field + Z" satisfies Eq. (j^^ . In spite 
of this non-uniqueness, we specify the components of the 



tensor Hab as Eq. (|5.4p , which is the gauge- invariant part 
of the linear-order metric perturbation associated with 
the longitudinal gauge. 

The non-uniqueness of the definitions of gauge- 
invariant variables is related to the "gauge-fixing" for 
the linear-order metric perturbations. Due to this non- 
uniqueness, we can consider the gauge- fixing in the first- 
order metric perturbation from two different points of 
view. The first point of view is that the gauge-fixing 
is to specify the gauge- variant part X". For example, 
the longitudinal gauge is realized by the gauge fixing 
X'^ — 0. Due to this gauge fixing X° = 0, we can regard 
the fact that perturbative variables in the longitudinal 
gauge are the completely gauge fixed variables. On the 
other hand, we may also regard that the gauge fixing is 
the specification of the gauge-invariant vector field Z° in 
Eq. (|5.5p . In this point of view, we do not specify the 
vector field X'^. Instead, we have to specify the gauge- 
invariant vector Z" or equivalently to specify the gauge- 
invariant metric perturbation ICab without specifying X° 
so that the first-order metric perturbation hab coincides 
with the gauge-invariant variables ICab when we fix the 
gauge X" so that X" + Z" = 0. These two different 
point of view of "gauge fixing" is equivalent with each 
other due to the non-uniqueness of the definition (|5.5p of 
the gauge-invariant variables. 



B. First-order Einstein equations 

Here, we derive the linear-order Einstein equation 
p.32p . To derive the components of the gauge invariant 
part of the linearized Einstein tensor '^^^Ga'' [H], which is 
defined by Eqs. (|3.8p . we first derive the components of 
the tensor H^^i^ [H] , which is defined in Eq. (|3.1ip with 
Aab — Hab and its component (j5.4p . These components 
are summarized in Ref. (Tl| . 

From Eq. (|3.8p . the component of '•^•'5^'' [H] are sum- 
marized as 



(1)0," [H] 
^'^G.'iH] 

^'^G,' [H] 



If. X , 

-GHdr, + 2A + 6K) * -6H^ $ 



2dr,D, * +2HD, <i> --{A + 2K) V, 



1 f (1) (1) 1 "^^^ 

— <^ 2d„D' m +2HD' $ +- (-A + 2K + AH^ - Ad,^H) 



(5.6) 
(5.7) 

(5.8) 



1 



D^D^ {% - $ ^ + I (-A + 2dl + 4Ha„ - 2K) % + {2Hd^ -I- Ad^H 2H^ A) $ | 7, 



(1) 



{dl^2Hdr,^2K- A) X, 



(5.9) 
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Straightforward calculations show that these components 
of the first-order gauge invariant perturbation '^^^G^'^ ['H] 
of the Einstein tensor satisfies the identity p.lSp . Al- 
though this confirmation is also possible without specifi- 
cation of the tensor Hab, the confirmation of Eq. (I3.13P 
through the explicit components (|5.6p - (|5.9p implies that 
we have derived the components of ^^•'5^'' [H] consis- 



tently. 

Next, we summarize the first-order perturbation of the 
energy momentum tensor for a scalar field. Since, at the 
background level, we assume that the scalar field (p is 
homogeneous as Eq. (|4.4p . the components of the gauge- 
invariant part of the first-order energy-momentum tensor 
'-^'Ta'' are given by 



(i>r ri 

' 71 



(Ih- 



1 / . 

— DVi -I- {d^(p) 



a —ipi 
dip 



(5.10) 
(5.11) 



r 



The second equation in (|5.1ip shows that there is no 
anisotropic stress in the energy-momentum tensor of the 
single scalar field. Then, we obtain 



Eq. (|4.8I) . We also note that only two of these equa- 
tions are independent. Further, the vector part of the 
component (^'^^ [H] = 87rG(i>7; shows that 



(1) (1) 
$ = * . 



(5.12) 



From Eqs. (|5.6p - (|5.1ip and (|5.12p . the components of 
scalar parts of the linearized Einstein equation p.32p are 
given as0 



(1) 



(A - 3Ha,, + 4i^ - dr^n - 2n^) $ 

= 4:7rG ( d„ipid„ip + a?'^ipi 
V dip 



(1) (1) 



(1) 



dV 



(5.13) 
(5.14) 



= AttG d„ipid„(p ~ a ——tpi . (5.15) 
V d'P J 

In the derivation of Eqs. (|5.13p - (|5.15p . we have used 



(1) 



(5.16) 



(1) 

The equation for the tensor mode Xij is given by 



(1) 



{d^ + 2nd,j + 2K - A) x/ = 0. 



(5.17) 



Combining Eqs. (|5.13p and (|5.15l) . we eliminate the 
potential term of the scalar field and thereby obtain 



(1) 



(a^ + A + AK) $ = 8T:Gdrjipidrjip 



(5.18) 



Further, using Eq. (|5.14p to express d,jipi in terms of 
(1) (1) 

d,j $ and <I) , we also eliminate 9,j</?i in Eq. (j5.18p . 
Hence, we have 



2\n- 



2_dy 



A-AK + 2 




(5.19) 



This is the master equation for the scalar mode pertur- 
bation of the cosmological perturbation in universe filled 
with a single scalar field. It is also known that Eq. (I5.19P 
reduces to a simple equation through a change of vari- 
ables 0. 



C. First-order Klein-Gordon equations 

Next, we consider the first-order perturbation of the 
Klein-Gordon equation (|3.28l) . By the straightforward 
calculations using Eqs. (HTTI) . ((5^ . (li^ . and the 

components -ff^j"*^ summarized in Ref. JTj, the gauge- 
(1) 

invariant part C(^k) of the first-order Klein-Gordon equa- 
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tion defined by Eq. (|3.28l) is given by 

/ (1) 
- ( 9,, $ +39,, \If I 

^ 2 X 2 5^1^. ^ 

+2a2 $ _(<^)+aVi^M 

= 0. (5.20) 

Througli tlie background Einstein equations (|4.6I) . 
(|47)) . and the first-order perturbations (|5.14p and (j5.19p 
of the Einstein equation, we can easily derive the 
first-order perturbation of the Klein-Gordon equation 
(|5.20[) 13]. Hence, the first-order perturbation of the 
Klein-Gordon equation is not independent of the back- 
ground and the first-order perturbation of the Einstein 
equation. Therefore, from the viewpoint of the Cauchy 
problem, any information obtained from the first-order 
perturbation of the Klein-Gordon equation should also 
be obtained from the set of the background and the first- 
order the Einstein equation, in principle. 

VI. EQUATIONS FOR THE SECOND-ORDER 
COSMOLOGICAL PERTURBATIONS 

Now, we develop the second-order perturbation the- 
ory on the cosmological background spacetime in Sec. ITVl 
within the general framework of the gauge-invariant per- 
turbation theory reviewed in Sec. |TT1 Since we have al- 
ready confirm the important step of our general frame- 
work, i.e., the assumption for the decomposition (I2.23P 
of the linear-order metric perturbation is correct. Hence, 
the general framework reviewed in Sec. is applica- 
ble. Applying this framework, we define the second-order 
gauge invariant variables of the metric perturbation in 
Sec. I VI Al In Sec. IVIBl we summarize the explicit com- 
ponents of the gauge invariant parts of the second-order 
perturbation of the Einstein tensor. In Sec. IVI C[ we 
summarize the explicit components of the second-order 
perturbation of the energy-momentum tensor and the 
Klein-Gordon equations. Then, in Sec. IVI Dl we derive 
the second-order Einstein equations in terms of gauge- 
invariant variables. The resulting equations have the 
source terms which constitute of the quadratic terms of 
the linear-order perturbations. Although these source 
terms have complicated forms, we give identities which 
comes from the consistency of all the second-order per- 
turbations of the Einstein equation and the Klein-Gordon 
equation in Sec. IVI El 

A. Gauge-invariant metric perturbations 

First, we consider the components of the gauge invari- 
ant variables for the metric perturbation of second order. 
The variable Lab defined by Eq. (|2.27p is transformed as 



Eq. (|2.28l) under the gauge transformation and we may 
regard the generator (Ja defined by Eq. (|2.29p as an ar- 
bitrary vector field on TWq from the fact that the gen- 
erator in Eq. (|2.29p is arbitrary. We can apply the 
procedure to find gauge invariant variables for the first- 
order metric perturbations (|5.4p in Sec. IV Al Then, we 
can accomplish the decomposition (|2.30l) . Following to 
the same argument as in the linear case, we may choose 
the components of the gauge invariant variables Cab in 
Eq. (jO^ as 

Cab = -2a^^<i {d'q)a{d'q)b + 2a^^Vi {dri)(^a{dx'')b) 

W (^-2 % 7,,+ % ^ {dx')a{dx^)b, (6.1) 

(2) (2) 

where v i and Xij satisfy the equations 

■ (2) ■ (2) 

(2) (2) 

The gauge invariant variables $ and vl> are the scalar 

(2) (2) 

mode perturbations of second order, and Vi and Xij ^re 
the second-order vector and tensor modes of the metric 
perturbations, respectively. 

Here, we also note the fact that the decomposition 
(|2.32p is not unique. This situation is similar to the case 
of the linear-order, but more complicated. In the defini- 
tion of the gauge invariant variables of the second-order 
metric perturbation, we may replace 

X'' = X''' - Z'", (6.3) 

where Z " is gauge invariant and A" " is transformed as 

yX''^ - xX'^ = (6.4) 

under the gauge transformation X\ ~^y\. This Z ° may 
be different from the vector Z° in Eq. (|5.5p . By the 
replacement (|6.3p . the second-order metric perturbation 
(|2.32p is given in the form 

^ab ='■ Jab + '2£x'hab + {£y' — -^X') 9ab, (6.5) 

where we defined 

Jab '■— Cab — £w9ab ~ 2£z'ICab 

-2£z'£zgab + £%'gab, (6.6) 

Y'" := Y" + + [X',ZT- (6.7) 

Here, the vector field W"" in Eq. (|6.7p constitute of some 
components of gauge invariant second-order metric per- 
turbation Cab like Z° in Eq. (|5.5p . The tensor field Jab 
is manifestly gauge invariant. The ^auge transformation 
rule of the new gauge- variant part y ° of the second-order 
metric perturbation is given by 

yY'^-xY'^ = e"2) + K(i),^r- (6.8) 
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Although Eq. (|631) is similar to Eq. (I02)) . the ten- 
sor fields Cab and Jab are different from each other. 
Thus, the definition of the gauge invariant variables for 
the second-order metric perturbation is not unique in 
a more complicated way than the linear order. This 
non-uniqueness of gauge-invariant variables for the met- 
ric perturbations propagates to the definition (j2.34p and 
(I2.35P of the gauge invariant variables for matter fields. 

In spite of the existence of infinitely many definitions 
of the gauge invariant variables, in this paper, we con- 
sider the components of Cab given by Eq. (|6.ip . Eq. (|6.ip 
corresponds to the second-order extension of the longitu- 
dinal gauge, which is called Poisson gauge = Y"' = 0. 



B. Einstein tensor 

Here, we evaluate the second-order perturbation of the 
Einstein tensor p. 71) with the cosmological background 



(OD . We evaluate the term (i)^/ [C] and (2)^^'' [H^U] 
in the Einstein equation p.33p . 

First, we evaluate the term ^^^Qa^ [C] in the Einstein 
equation p. 331) . Because the components (|6.ip of Cab 
are obtained through the replacements 



(1) (2) (1) (2) 
$ ^ $ , Vi ^ V 



(1) (2) 



x'ij ^ Xy- (6-9) 



in the components (|5.4p of T-Lab, we easily obtain the com- 
ponents of ^^^Qj'lC] through the replacements (|6.9p in 

Eqs. (ESD-dSH). 

From Eq. (|5.4I) . we can derive the components of 
i^)gj> = {2)g^b^n,n] defined by Eqs. ^^ISJJ^ in 
a straightforward manner. Here, we use the results 
(|5.12p and (|5.16p of the first-order Einstein equations, 
for simplicity. Then the explicit components '•^-'^/a'' = 
^^'5a''[H,H] are summarized as 



(2)C V 



(2)^ 



2 r (1) ,.(1) (1) , (1) , n ^ (l)^Vfc 

— -3i?fc $ i?^ $ -8 $ A $ -3 fa^ $ j -12(H2 + X)f$j +A£'fe $x 

1 ni 1 iw 1 i\\ 1 i-w 



(1) 
$ 



(1) (1) (1) (1) / (1) (1) 1 (1) /i\ >■ 



(1) 



(1) (1) (1) (1) ■ (1) (1) 1 ^t- (1) ^u- (1) 

A $ -2D^ $ 9^ $ -t-D-' $ 9^ Xij -driD^ $Xu -^^r, X A Xfcj + X X{\k 



Xkl drj 

W (1) 



.k]l 



(6.10) 
,(6.11) 
. (6.12) 



"f (1) ,(1) (1), (1) (1) (1) (1) (1) 

<^ -iDk $ D'^ $ -4 $ (A + if ) $ -d^ $ $ -87^ $ $ -4 (29^^ + H 

(1) , (1) (1) , (1) , , , (1) (1) 

+2D., $ LIJ $ +4 $ D,D^ <^ +x^' {dl + 2ndr,) $ +Dk $ 



/'(l) 
2^ ' $ 



>7» 



X^' Xik 



(1) 



(1) 



(1) (1) (1) ^'1. ■ (1) (1) (1) 

-2D'' $ Dk X,-' -2 $ (A - 2if) X,' -A $x.' ^x''' $x™ -DiDk $x 7, 

S a, x'^^' S D^' x'l^' +^A' x'i A x"= +^ AA X™' 



(1) 



(1) 



(1) 



(i)(i) 



1 



(1) 



(1) I (1) 



1 (1) n) 1 (1). 1 ri) 

-- x"" i^iA +- x'" D„M X,' ~2 + ^^^"^ - A + 2X) 



2 

1 j 3 (1) 'fc; , , 

'214 '' ^ + ^'^^ 

hA ;^:ii?['x^'^l'47 



(1) 



{ii {d^ + 2Hd, -A + K) x'"" ~^Dk xiL D'^ x™' 



(6.13) 



We have checked the identity p.l4p through Eqs. (|6.10p - 



(|6.13p . Then, we may say that the expressions (|6.10p - 
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(|6.13|) are self-consistent. 



C. Energy-momentum tensor and Klein-Gordon 
equation 

Here, we summarize the explicit components of the 
gauge-invariant parts (j3.23p of the second-order pertur- 



bation of energy momentum tensor for a single scalar 
field in terms of gauge-invariant variables. Through 
Eqs. (j44|) . (|6TT|) . the components of Eq. (|X23l) are 

derived by the straightforward calculations. In this pa- 
per, we just summarize the components of in the 
situation where the first-order Einstein equations ()5.12p 
and (|5.16p are satisfied: 



(1) 



(1) . 



(1) 



-DiifiD ifii-a ^2-Q^ - a [tpi) -^-j 



(2) 



(6.14) 

(6.15) 
(6.16) 

(6.17) 



r 



More generic formulae for the components of '•^-'7^'' are 
given in Ref. [l^ . 

Next, we show the gauge-invariant second-order the 
Klein-Gordon equation. We only consider the simple sit- 
uation where Eqs. (|5.12p and (|5.16p are satisfied. The 
formulae for more generic situation is given in Ref. 
Through Eqs. (lO)) . (|6T|) . (|44|) . the second-order per- 
turbation of the Klein-Gordon equation (I3.29p is given 

by 



(2) 



= df^^2 + 2nd^^2-A^2 

(2) (2)\ 



-2a^ $ — (^) +a>2 
dip 



(^) 



0, 



(6.18) 



where we defined 
^ (1 



$ d,,i^i -I- 8 <J) Aipi - 4a $ <y5i-^r:7('^) 



d^V (1) (1) 

-a ((pi) -^:r^{f) + 8 ^ ^ dr^(p 



(1) 



(1) 



(1) 



-2 Avi + dr,^ x'' xi, ■ (6.19) 

In Eq. (|6.18p . "^{K) is the source term which is the 
collection of the quadratic terms of the linear-order 



perturbations in the second-order perturbation of the 
Klein-Gordon equation. If we ignore this source term, 
Eq. (j6.18l) coincide with the first-order perturbation of 
the Klein-Gordon equation. From this source term ()6.19p 
of the Klein-Gordon equation, we can see that the mode- 
mode coupling due to the non-linear effects appear in the 
second-order Klein-Gordon equation. 

We cannot discuss solutions to Eq. (|6.18p only through 
this equation, since this includes metric perturbations. 
To determine the behavior of the metric perturbations, 
we have to treat the Einstein equations simultaneously. 
The second-order Einstein equation is shown in Sec lVIP] 



D. Einstein equations 

Here, we show the all components of the second-order 
Einstein equation (I3.33p . All components of Eq. p.33p 
are summarized as 



X (2) 



dV 



, (2) 

K) $ 



-47rG ydn^dr^V2 + a V^-q^J = Tq, 

(2) (2) 1 (2) 

2dnD, 4- +2HA $ - - (A + 2X) W 
~87rGDi(p2dn(p = E^, 



(6.20) 



(6.21) 
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/(2) {2)\ { (2) , ^ „ ^ , . N (2)1 



dV_ 

dtp 



^49.) ( t^i ) + \{dl + 2Hd^ + 2K-A) x\j -SnG (d,j^d^(p2 - a^(p2^{^) ) = T.^ , (6.22) 



where Fq, Ti Tij are the collection of the quadratic term of the first-order perturbations as follows: 
To := 4^Gf (a^(^i)2+A¥'i£'Vi+a2((^i)2_j - 4d^n (<P j - 2 $ ^2 $ -3A $ i?^ $ -10 $ A $ 

-3 i^d, $y - 16X - sn' + A A $x'^ +ld, d, x^' +H d, x'^ 

3 (1) 1 fi) 1 '■^^ (1) 1 (1) 

»™ i?'^ X"' Xi!n x"'^ ^2 x'" A XiL iliLx'"'- (6.23) 

(1) (1) (1) (1) (1) (1) ■ (1) (1) ■ (1) (1) 

r, := l&TTGdrj^iDiipi -45^ $ A $ +8H $ A $ -8 $ 9^ A $ +2 A $ 3^ x/« -23,, A $x4 

- 2 5, x'^ - Xfc; d^D, x'' + 9,^?^ 'xL ; (6.24) 
r,, := 16^GA</5iA</5i + 87rG |(a,,^i)2 - Di^iD^^i - a2((^i)2|Z| _ 4^)^ A "8 AA 

/ (1) (1) (1) (1) ^ 9/(l)\^ (1)^ (1) ^9 

+ ezjfc $ z?*^ $ +4 $ A $ +2 f a,, $ J + 8dr,n ( * ) + i6-h2 f $ J + 16H $ 9,, $ -4 $ 92 $ 

„ „9(1)(1) t (1) t (1) . (1) t (1) 

-4-H9,, $X»j -29' $X»j -4i?' $ D(^ Xj)k +4i?' $ Xrj ~&K $x,j +4 $ A x4 -4Z)^i?(, $x,)fc 

(1) (1) (1) (1) n~) ("^^ fi) (1) 1 (I'l 

+2A $^4- +2 A A <i>x"^ 7., + x.fe 9,, x/ -D" Xu A x/ xL D' Xjk -^D, x"= A Xifc 

- X/i A A X'"' +2 x'" ^/^(. X,ln - x'" i?™^; X^i 

-4 I 39, xii drj x" -3 A XiL X™' +2 A XiL x'"'^ -4if x/ix"" 1 1.,- (6.25) 

Here, we used Eqs. gill), (lET^ . (Ol) . (lETCl) and (EH]). 
The tensor part of Eq. (|6.22p is given by 

(92 + 2-W9,, + 2X - A) % - 2r., - ^7,,r/^ - 3 (^AA - (A + {^/^-'D^DiVf! - ip/'^ 

+4 {i:'(,(A + 2K)-^Dj-jA-^D^DkTi - D^,{A + 2K)-^D''rj-jk} ■ (6.26) 

This tensor mode is also called the second-order gravitational waves. 

Further, the vector part of Eqs. (|6.2ip and (|6.22p yields the initial value constraint and the evolution equation of 

(2) 

the vector mode Vj : 



A + 2K 



- { A A-i AEfe - r,} , 9, (^a2 ^^l^ = { AA-i A AEfc' - AE, ^} . (6.27) 
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Finally, scalar part of Eqs. (I6.21[) - (j6.22p are summarized as 

(2) (2) 



2d^ +2n $ -87rG^25,,^ = A-'OTk, (6.28) 

*| - ^(A + f^A-'DW,r.^^ - ir,'^! , (6.29) 

-d^ - 5Hd„ + + %^ - (2d„H + Hd,^ + + ^A^ $ _87rG'aV2^ = To - (6.30) 

= -ro - ^r,' + A-^D'DjTj + id, - A-'D'r, 



where T^^ := j'^^Tik and T^^'^ = -f'^T.j. Eq. (lOTIl is the 
second-order extension of Eq. (I5.19p . which is the master 
equation of scalar mode of the second-order cosmological 
perturbation in a universe filled with a single scalar field. 

Thus, we have a set of ten equations for the second- 
order perturbations of a universe filled with a single scalar 
field, Eqs. (|6.26P " (|6.3ip . To solve this system of equa- 
tions of the second-order Einstein equation, first of all, 
we have to solve the linear-order system. This is ac- 
complished by solving Eq. (15. 19^ to obtain the potential 
(1) (1) 
$ , (^1 is given through (|5.14p . and the tensor mode X ij 

is given by solving Eq. (|5.17p . Next, we evaluate the 
quadratic terms, Fq, F^ and Tij of the linear-order pertur- 
bations, which are defined by Eqs. (|6.23p - (|6.25l) . Then, 
using the information of Eqs. (I6.23l) - (|6.25p . we estimate 
the source term in Eq. (j6.3ip . If we know the two in- 
dependent solutions to the linear-order master equation 
(|5.19p . we can solve Eq. (|6.3ip through the method us- 
ing the Green functions. After constructing the solution 

(2) 

$ to Eq. (|6.3ip . we can obtain the second-order met- 

(2) 

ric perturbation vj/ through Eq. (|6.29l) . Thus, we have 
obtained the second-order gauge invariant perturbation 
(f2 of the scalar field through Eq. (|6.28p . Thus, the all 

(2) (2) 

scalar modes $ , v]/ , (^2 are obtained. Equation (|6.30p 
is then used to check the consistency of the second-order 
perturbation of the Klein Gordon equation (|6.18p as in 
Sec.lVm 

(1) 

For the vector-mode, Vi of the first-order identically 
vanishes due to the momentum constraint (j5.16p for the 
linear-order metric perturbations. On the other hand, 
in the second-order, we have evolution equation (j6.27p 

(2) . 

of the vector mode Ui with the initial value constraint. 
This evolution equation of the second-order vector mode 
should be consistent with the initial value constraint, 
which is confirmed in Sec. rvTEl Equations (pTTf)) also im- 



ply that the second-order vector-mode perturbation may 
be generated by the mode couplings of the linear order 
perturbations. As the simple situations, the generation 
of the second-order vector mode due to the scalar-scalar 
mode coupling is discussed in Refs. f26j . 

The second-order tensor mode is also generated by the 
mode-coupling of the linear-order perturbations through 
the source term in Eq. (|6.26p . Note that Eq. (|6.26p is al- 
most same as Eq. (|5.17p for the linear-order tensor mode, 
except for the existence of the source term in Eq. (|6.26l) . 
If we know the solution to the linear-order Einstein equa- 
tions (|5.17p and (|5.19p . we can evaluate the source term 
in Eq. (|6.26p . Further, we can solve Eq. (|6.26p through 
the Green function method. This leads the generation of 
the gravitational wave of the second order. Actually, in 
the simple situation where the first-order tensor mode ne- 
glected, the generation of the second-order gravitational 
waves discussed in some literature [27|]. 

E. Consistency of equations for second-order 
perturbations 

Now, we consider the consistency of the second-order 
perturbations of the Einstein equations (|6.28p - (|6.3ip for 
the scalar modes, Eqs. (|6.27p for vector mode, and the 
Klein- Gordon equation (16.181) . The consistency check of 
these equations are necessary to guarantee that the de- 
rived equations are correct, since the second-order Ein- 
stein equations have complicated forms owing to the 
quadratic terms of the linear-order perturbations that 
arise from the nonlinear effects of the Einstein equations. 

Since the first equation in Eqs. (j6.27p is the initial value 

(2) 

constraint for the vector mode and it should be consis- 
tent with the evolution equation, i.e., the second equation 
of Eqs. (|6.27p . these equations should be consistent with 
each other from the general arguments of the Einstein 
equation. Explicitly, these equations are consistent with 
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each other if the equation 

dr^Tk + 2HVk - D'Tik = (6.32) 

is satisfied. Actually, through the first-order perturbative 
Einstein equations (|5.f 41) . (I5.f 91) . (I5.f 7p . we can confirm 
the equation (|6.32p . This is a trivial result from a general 
viewpoint, because the Einstein equation is the first class 
constrained system. However, this trivial result implies 
that we have derived the source terms Ti and Tij of the 
second-order Einstein equations consistently. 

Next, we consider Eq. (|6.30p . Through the second- 
order Einstein equations (|08l) . ([09| . (|63T|) . and the 
background Klein-Gordon equation (|4.9|) . we can con- 
firm that Eq. (|6.30p is consistent with the set of the 
background, first-order and other second-order Einstein 
equation if the equation 

(a^ + 2%) D^Tk - DWT.j = (6.33) 

is satisfied under the background and first-order Einstein 
equations. Actually, we have already seen that Eq. (I6.32p 
is satisfied under the background and first-order Einstein 
equations. Taking the divergence of Eq. (|6.32|) . we can 
immediately confirm Eq. (|6.33p . Then, Eq. (|6.30p gives 
no information. 

Thus, we have seen that the derived Einstein equations 
of the second order (|6.27l) - (|6.3ip are consistent with each 
other through Eq. (|6.32p . This fact implies that the de- 
rived source terms Ti and of the second-order per- 
turbations of the Einstein equations, which are defined 
by Eqs. (|6.24p and (|6.25l) . are correct source terms of the 
second-order Einstein equations. On the other hand, for 
To, we have to consider the consistency between the per- 
turbative Einstein equations and the perturbative Klein- 
Gordon equation as seen below. 

Now, we consider the consistency of the second-order 
perturbation of the Klein-Gordon equation and the Ein- 
stein equations. The second-order perturbation of the 
Klein-Gordon equation is given by Eq. (I6.18P with the 

(2) 

source term (|6.19p . Since the vector mode Vi and tensor 

(2) 

mode X ij of the second-order do not appear in the ex- 
pressions (|6.18p of the second-order perturbation of the 
Klein-Gordon equation, we may concentrate on the Ein- 
stein equations for scalar mode of the second order, i.e., 
Eqs. (11211), and (|OT|) with the definitions ((l!^ - 

(|6.25p of the source terms. As in the linear case, the 
second-order perturbation of the Klein-Gordon equation 
should also be derived from the set of equations consisting 
of the second-order perturbations of the Einstein equa- 
tions ((Ogi) . (jOTI) . the first-order perturbations 
of the Einstein equations (|5.12p . (|5.14p . (15.191) . and the 
background Einstein equations (14.61) and (|4.7p . Actually, 
from these equation, we can show that the second-order 
perturbation of the Klein-Gordon equation is consistent 
with the background and the second-order Einstein equa- 
tions if the equation 

2 {dr, +n)TQ- D'^Tk + UT^'^ 

+8T:GdnipE(^K) = (6.34) 



is satisfied under the background and the first-order Ein- 
stein equations. Further, we can also confirm Eq. (|6.34p 
through the background Einstein equations (|4.6p and 
(|4.7p . the scalar part of the first-order perturbation of 
the momentum constraint (|5.14p . the evolution equations 
(I5.19P and (I5.17P for the first order scalar and tensor 
modes in the Einstein equation. 

As shown in Ref. ,13 ], the first-order perturbation of 
the Klein-Gordon equation is derived from the back- 
ground and the first-order perturbations of the Einstein 
equation. In the case of the second-order perturbation, 
the Klein-Gordon equation (|6.18p can be also derived 
from the background, the first-order, and the second- 
order Einstein equations. The second-order perturba- 
tions of the Einstein equation and the Klein- Gordon 
equation include the source terms Fq, F^, F^j, and '^{k) 
due to the mode-coupling of the linear-order pertur- 
bations. The second-order perturbation of the Klein- 
Gordon equation gives the relation (j6.34p between the 
source terms Fq, F^, F^ , "^(k) and we have also confirmed 
that Eq. (|6.34p is satisfied due to the background, the 
first-order perturbation of the Einstein equations, and 
the Klein-Gordon equation. Thus, the second-order per- 
turbation of the Klein- Gordon equation is not indepen- 
dent of the set of the background, the first-order, and the 
second-order Einstein equations if we impose on the Ein- 
stein equation at any conformal time rj. This also implies 
that the derived formulae of the source terms Fq, F^, F^ , 
and are consistent with each other. In this sense, 

we may say that the formulae (j6.23p - (j6.25p and (|6.19p 
for these source terms are correct. 



VII. SUMMARY AND DISCUSSIONS 

In this paper, we summarized the current status of the 
formulation of the gauge-invariant second-order cosmo- 
logical perturbations. Although the presentation in this 
paper is restricted to the case of the universe filled by a 
single scalar field, the essence of the general framework 
of the gauge-invariant perturbation theory is transpar- 
ent through this simple case. The general framework 
of the general relativistic higher-order gauge-invariant 
perturbation theory can be separated into three parts. 
First one is the general formulation to derive the gauge- 
transformation rules (|2.18p and (|2.19p . Second one is the 
construction of the gauge-invariant variables for the per- 
turbations on the generic background spacetime inspect- 
ing gauge-transformation rules (I2.18P and (|2.19p and the 
decomposition formula (|2.36l) and (|2.37l) for perturba- 
tions of any tensor field. Third one is the application 
of the above general framework of the gauge-invariant 
perturbation theory to the cosmological situations. 

To derive the gauge-transformation rules (|2.18p and 
(|2.19p . we considered the general arguments on the Tay- 
lor expansion of an arbitrary tensor field on a manifold, 
the general class of the diffeomorphism which is wider 
than the usual exponential map, and the general for- 
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mulation of the perturbation theory. This general class 
of diffconiorphism is represented in terms of the Tay- 
lor expansion (j2.2p of its pull-back. As commented in 
Sec. Ill Al this general class of diffeomorphism does not 
form a one-parameter group of diffeomorphism as shown 
through Eq. (|2.3p . However, the properties p.3|) do not 
directly mean that this general class of diffeomorphism 
does not form a group. One of the key points of the prop- 
erties of this diffeomorphism is the non-commutativity of 
generators f " and ^2 of each order. Although the expres- 
sion of the n-th order Taylor expansion of the pull-back of 
this general class is discussed in Ref. [22|, when we con- 
sider the situation of the n-th order perturbation, this 
non-commutativity becomes important 0. Therefore, to 
clarify the properties of this general class of diffeomor- 
phism, we have to take care of this non-commutativity of 
generators. Thus, there is a room to clarify the properties 
of this general class of diffeomorphism. 

Further, in Sec. Ill CI we introduced a gauge choice 
X\ as an exponential map, for simplicity. On the other 
hand, we have the concept of the general class of dif- 
feomorphism which is wider than the class of the ex- 
ponential map. Therefore, we may introduce a gauge 
choice as one of the element of this general class of dif- 
feomorphism. However, the gauge-transformation rules 
(|2.18[) and (|2.19p will not be changed even if we general- 
ize the definition of a each gauge choice as emphasized in 
Sec. Ill Cl Although there is a room to sophisticate in log- 
ical arguments to derive the gauge-transformation rules 
(|2.18[) and (|2.19p , these are harmless to the development 
of the general framework of the gauge-invariant pertur- 
bation theory shown in Sees. Ill C| III El IIIH and their 
application to cosmological perturbations in Sec. IIVI 

As emphasize in Sec. Ill E| our starting point to con- 
struct gauge invariant variables is the assumption that 
we already know the procedure for finding gauge in- 
variant variables for the linear metric perturbations as 
Eq. i2. 23\) . This is highly nontrivial assumption on a 
generic background spacetime. The procedure to accom- 
plish the decomposition (|2.23p completely depends on 
the details of the background spacetime. In spite of this 
non-triviality, this ass ump tion is almost correct in some 
background spacetime [28|. Further, once we accept this 
assumption, we can develop the higher-order perturba- 
tion theory in an independent manner of the details of 
the background spacetime. We also expect that this gen- 
eral framework of the gauge-invariant perturbation the- 
ory is extensible to n-th order perturbation theory, since 
our procedure to construct gauge-invariant variables can 
be extended to the third-order perturbation theory with 
two-parameter 0. Due to this situation, in Ref. [l3|, we 
propose the conjecture which state that the above as- 
sumption for the decomposition of the linear-order met- 
ric perturbation is correct for any background spacetime. 
We may also say that the most nontrivial part of our gen- 
eral framework of higher-order gauge-invariant perturba- 
tion theory is in the above assumption. Further, as em- 
phasized in Sec. IV Al we assumed the existence of some 



Green functions to accomplish the decomposition (j2.23p 
and this assumption exclude some perturbative modes of 
the metric perturbations from our consideration, even in 
the case of cosmological perturbations. For example, ho- 
mogeneous modes of perturbations are excluded in our 
current arguments of the cosmological perturbation the- 
ory. These homogeneous modes will be necessary to dis- 
cuss the comparison with the arguments based on the 
long-wavelength approximation. Therefore, we have to 
say that there is a room to clarify even in the cosmolog- 
ical perturbation theory. 

Even if the assumption is correct on any background 
spacetime, the other problem is in the interpretations 
of the gauge-invariant variables. We have commented 
on the non-uniqueness in the definitions of the gauge- 
invariant variables through Eqs. (|5.5p and (|6.5p . This 
non-uniqueness in the definition of gauge-invariant vari- 
ables also leads some ambiguities in the interpretations of 
gauge-invariant variables. On the other hand, as empha- 
size in Sec. Ill CI any observations and experiments are 
carried out only on the physical spacetime through the 
physical processes on the physical spacetime. For this 
reason, any direct observables in any observations and 
experiments should be independent of the gauge choice. 
Further, the non-uniqueness in the definitions the gauge- 
invariant variables expressed by Eqs. (|5.5p and (|6.5p have 
the same form as the decomposition formulae (|2.36p and 
(I2.37p . Therefore, if the statement that any direct ob- 
servables in any observations and experiments is indepen- 
dent of the gauge choice is really true, this also confirm 
that the non-uniqueness of the definition of the gauge- 
invariant variables also have nothing to do with the direct 
observables in observations and experiments. These will 
be confirmed by the clarification of the relations between 
gauge-invariant variables and observables in experiments 
and observations. To accomplish this, we have to spec- 
ify the concrete process of experiments and observations 
and clarify the problem what are the direct observables 
in the experiments and observations and derive the rela- 
tions between the gauge-invariant variables and observ- 
ables in concrete observations and experiments. If these 
arguments are completed, we will be able to show that 
the gauge degree of freedom is just unphysical degree of 
freedom and the non-uniqueness of the gauge-invariant 
variables have nothing to do with the direct observables 
in the concrete observation or experiment, and then, we 
will be able to clarify the precise physical interpretation 
of the gauge-invariant variables. 

For example, in the case of the CMB physics, we 
can easily see that the first-order perturbation of the 
CMB temperature is automatically gauge invariant from 
Eq. (I2.36p . because the background temperature of CMB 
is homogeneous. On the other hand, the decomposi- 
tion formula (|2.37p of the second order yields that the 
theoretical prediction of the second-order perturbation 
of the CMB temperature may depend on gauge choice, 
since we do know the existence of the first-order fluctu- 
ations as the temperature anisotropy in CMB. However, 
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as emphasized above, the direct observables in observa- 
tions should be gauge invariant and the gauge-variant 
term in Eq. (j2.37|) should be disappear in the direct ob- 
servables. Therefore, we have to clarify the how gauge- 
invariant variables are related to the observed tempera- 
ture fluctuations and the gauge-variant term disappear 
in the observable. 

Although there are some rooms to accomplish the com- 
plete formulation of the second-order cosmological per- 
turbation theory, we derived all the components of the 
second-order perturbation of the Einstein equation with- 
out ignoring any types modes (scalar-, vector-, tensor- 
types) of perturbations in the case of a scalar field system. 
In our formulation, any gauge fixing is not necessary and 
we can obtain all equations in the gauge-invariant form, 
which are equivalent to the complete gauge fixing. In 
other words, our formulation gives complete gauge-fixed 
equations without any gauge fixing. Therefore, equations 
obtained in a gauge-invariant manner cannot be reduced 
without physical restrictions any more. In this sense, the 
equations shown here are irreducible. This is one of the 
advantages of the gauge-invariant perturbation theory. 

The resulting Einstein equations of the second order 
show that any type of mode-coupling appears as the 
quadratic terms of the linear-order perturbations owing 
to the nonlinear effect of the Einstein equations, in prin- 
ciple. Perturbations in cosmological situations are clas- 
sified into three types: scalar, vector, and tensor. In 
the second-order perturbations, we also have these three 
types of perturbations as in the case of the first-order 
perturbations. Furthermore, in the equations for the 
second-order perturbations, there are many quadratic 
terms of linear-order perturbations owing to the non- 
linear effects of the system. Owing to these nonlinear 
effects, the above three types of perturbations couple 
with each other. In the scalar field system shown in this 
paper, the first-order vector mode does not appear due 
to the momentum constraint of the first-order perturba- 
tion of the Einstein equation. Therefore, we have seen 
that three types of mode-coupling appear in the second- 
order Einstein equations, i.e., scalar-scalar, scalar-tensor, 
and tensor-tensor type of mode coupling. In general, all 
types of mode-coupling may appear in the second-order 
Einstein equations. Actually, in Ref. [ll], we also de- 
rived the all components of the Einstein equations for a 
perfect fluid system and we can see all types of mode- 
coupling, i.e., scalar-scalar, scalar- vector, scalar-tensor, 
vector- vector, vector-tensor, tensor-tensor types mode- 
coupling, appear in the second-order Einstein equation, 
in general. Of course, in the some realistic situations of 
cosmology, we may neglect some modes. In this case, we 
may neglect some mode-coupling. However, even in this 
case, we should keep in mind the fact that all types of 
mode-couplings may appear in principle when we discuss 
the realistic situations of cosmology. We cannot deny 
the possibility that the mode-couplings of any type pro- 
duces observable effects when the quite high accuracy of 
observations is accomplished. 



Even in the case of the single scalar field discussed 
in this paper, the source terms of the second-order Ein- 
stein equation show the mode-coupling of scalar-scalar, 
scalar-tensor, and the tensor-tensor types as mentioned 
above. Since the tensor mode of the linear order is also 
generated due to quantum fluctuations during the infla- 
tionary phase, the mode-couplings of the scalar-tensor 
and tensor-tensor types may appear in the inflation. 
If these mode-couplings occur during the inflationary 
phase, these effects will depend on the scalar-tensor ratio 
r. If so, there is a possibility that the accurate observa- 
tions of the second-order effects in the fluctuations of the 
scalar type in our universe also restrict the scalar-tensor 
ratio r or give some consistency relations between the 
other observations such as the measurements of the B- 
mode of the polarization of CMB. This will be a new 
effect that gives some information on the scalar-tensor 
ratio r. 

Furthermore, we have also checked the consistency be- 
tween the second-order perturbations of the equations of 
motion of matter field and the Einstein equations. In 
the case of a scalar field, we checked the consistency 
between the second-order perturbations of the Klein- 
Gordon equation and the Einstein equations. Due to this 
consistency check, we have obtained the consistency re- 
lations between the source terms in these equations Fq, 
Ti, Tij, and '^{k)j which are given by Eqs. (|6.32l) and 
(I6.34p . We note that the relation (|6.32p comes from the 
consistency in the Einstein equations of the second order 
by itself, while the relation (|6.34p comes from the con- 
sistency between the second-order perturbation of the 
Klein-Gordon equation and the Einstein equation. We 
also showed that these relations between the source terms 
are satisfied through the background and the first-order 
perturbation of the Einstein equations in Ref. [isj . This 
implies that the set of all equations are self-consistent 
and the derived source terms Fq, F^, Tij, and '^{k) are 
correct. We also note that these relations are indepen- 
dent of the details of the potential of the scalar field. 

Thus, we have derived the self-consistent set of equa- 
tions of the second-order perturbation of the Einstein 
equations and the evolution equations of matter fields in 
terms of gauge-invariant variables. As the current status 
of the second-order gauge-invariant cosmological pertur- 
bation theory, we may say that the curvature terms in 
the second-order Einstein tensor p.33p . i.e., the second- 
order perturbations of the Einstein tensor, are almost 
completely derived although there remains the problem 
of homogeneous modes as mentioned above. After com- 
plete the problem of homogeneous modes, we have to 
clarify the physical behaviors of the second-order cosmo- 
logical perturbation in the single scalar field system in 
the context of the inflationary scenario. This is the pre- 
liminary step to clarify the quantum behaviors of second- 
order perturbations in the inflationary universe. Further, 
we also have to carry out the comparison with the re- 
sult by long-wavelength approximations. If these issues 
are completed, we may say that we have completely un- 
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derstood the properties of the second-order perturbation 
of the Einstein tensor. The next task is to clarify the 
nature of the second-order perturbation of the energy- 
momentum tensor through the extension to multi-fluid 
or multi-field systems. Further, we also have to extend 
our arguments to the Einstein Boltzmann system to dis- 
cuss CMB physics, since we have to treat photon and 
neutrinos through the Boltzmann distribution functions. 
This issue is also discussed in some literature 0, 13. If 
we accomplish these extension, we will be able to clarify 
the Non-linear effects in CMB physics. 

Finally, readers might think that the ingredients of this 
paper is too mathematical as Astronomy. However, we 
have to emphasize that a high degree of the theoretical 
sophistication leads unambiguous theoretical predictions 
in many case. As in the case of the linear-order cosmo- 
logical perturbation theory, the developments in observa- 
tions are also supported by the theoretical sophistication 
and the theoretical sophistication are accomplished mo- 
tivated by observations. In this sense, now, we have an 
opportunity to develop the general relativistic second- 
order perturbation theory to a high degree of sophistica- 
tion which is motivated by observations. We also expect 
that this theoretical sophistication will be also useful to 
discuss the theoretical predictions of Non-Gaussianity in 
CMB and comparison with observations. Therefore, I 
think that this opportunity is opened not only for obser- 
vational cosmologists but also for theoretical and math- 
ematical physicists. 
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Appendix A: Derivation of the generic 
representation of the Taylor expansion of tensors on 
a manifold 

In this section, we derive the representation of the co- 
efficients of the formal Taylor expansion (|2.2p of the pull- 
back of a diffeomorphism in terms of the suitable deriva- 
tive operators. The guide principle of our arguments is 
the following theorem [20l [23j . 

Theorem A.l. Let J) be a derivative operator acting on 
the set of all the tensor fields defined on a differentiate 
manifold M and satisfying the following conditions: (i) 
it is linear and satisfies the Leibniz rule; (ii) it is tensor- 
type preserving; (Hi) it commutes with every contraction 



of a tensor field; and (iv) it commutes with the exte- 
rior differentiation d. Then, T) is equivalent to the Lie 
derivative operator with respect to some vector field 
i.e., T) ~ £^. 

The prove of the assertion of Theorem lA.ll is given 
in Ref. [l^l as follows. When acting on functions, the 
derivative operator V defines a vector field ^ through the 
relation 

Vf^:^if)^£^f, VfeT{M). (Al) 

The assertion of the Theorem for an arbitrary tensor field 
is hold iff the assertions for an arbitrary scalar function 
and for an arbitrary vector field V are hold. To do this, 
we consider the scalar function V{f) and we obtain 

mif)) = aV{f)) (A2) 

through Eq. (jAII) . Through the conditions (i)-(iv) of 23, 
V{V{f)) is also given by 

V{V{f)) = V{df{V))^V{C{df®V)] 
= C{D{df®V)] 
= C {V{df) ®V + df® VV} 
= C{d{Vf)®V + df ®VV} 
= d{Vf){V)+df{VV) 
= V{Vf) + {VV){f) (A3) 

Then we obtain 

= i{v{f))^vm))^[£.,v]{f) 

= {£iV){f) (A4) 

for an arbitrary /, i.e., 

VV = £^V. (A5) 

Through Eqs. (|Aip and (|A5I) . we can recursively show 

VQ = £^Q (A6) 

for an arbitrary tensor field Q[23j. 

Now, we consider the derivation of the Taylor expan- 
sion (|2.ip . As in the main text, we first consider the 
representation of the Taylor expansion of for an ar- 
bitrary scalar function / G F{M): 

where J-{M) denotes the algebra of C°° functions on A4. 
Although the operator d/dX in the bracket {*}a=o of 
Eq. ()A7|) are simply symbolic notation, we stipulate the 
properties 
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for V/ G where n is an arbitrary finite integer. 

These properties imply that the operator d/dX is in fact 
not simply symbolic notation but indeed the usual par- 
tial differential operator on K. We note that the prop- 
erty (|A9p is the Leibniz rule, which plays important roles 
when we derive the representation of the Taylor expan- 
sion (jA7l) in terms of suitable Lie derivatives. 

Together with the property (IA9|) . Theorem I A . II yields 
that there exists a vector field so that 

Actually, the conditions (ii)-(iv) in Theorem I A . 1 1 are sat- 
isfied from the fact that $^ is the pull-back of a dif- 
feomorphism and (i) is satisfied due to the property 
(1X91) . 

Next, we consider the second-order term in Eq. (jA7p . 
Since we easily expect that the second-order term in 
Eq. (jA7l) may includes , we define the derivative op- 
erator £2 by 

{^('J'a/)}^_^ =■■ {q2)+a£l)f, (All) 

where a is determined so that C2 satisfy the conditions 
of Theorem lA.il The conditions (ii)-(iv) in Theorem lA.il 
for £2 are satisfied from the fact that $^ is the pull-back 
of a diffeomorphism $>. Further, £2 is obviously linear 
but we have to check £2 satisfy the Leibniz rule, i.e., 

£2 if) = 2/£2/ (A12) 

for V/ G T{M). To do this, we use the properties (lASp 
and (|A9|) . then we can easily see that the Leibniz rule 
(|Al2p is satisfied iff a = 1 and we may regard £2 as the 
Lie derivative with respect to some vector field. Then, 
when and only when a = 1, there exists a vector field ^2 
such that 

£2/ = £iJ (A13) 

and 

{^('J'a/)}^ „ + (A14) 

Thus, wc have seen that the Taylor expansion (jA7p for 
an arbitrary scalar function / is given by Eq. (|2.2p . 

Although the formula (12. 2p of the Taylor expansion is 
for an arbitrary scalar function, we can easily extend this 
formula to that for an arbitrary tensor field Q as the as- 
sertion of Theorem lA.il The proof of the extension of the 
formula (|2.2p to an arbitrary tensor field Q is completely 
parallel to the proof of the formula (|2.2I) for an arbitrary 
scalar function if we stipulate the properties 



instead of Eqs. (jASp and (jA9p . As the result, we obtain 
the representation of the Taylor expansion for an arbi- 
trary tensor field Q. 

Appendix B: Derivation of the perturbative Einstein 
tensors 

Following the outline of the calculations explained in 
Sec. IIII Al we first calculate the perturbative expansion 
of the inverse metric. The perturbative expansion of the 
inverse metric can be easily derived from Eq. (|2.22p and 
the definition of the inverse metric 

g'^'Sbc = S^. (Bi) 

We also expand the inverse metric 

gab 

the form 

-ab ^ gab ^ ^(l)-ab ^ ^X^-^Y^ (B2) 

Then, each term of the expansion of the inverse metric is 
given by 

(1)^'^'' = -h"^^ ^^^b = 2h'"'hJ' - (B3) 

To derive the formulae for the perturbative expansion 
of the Riemann curvature, we have to derive the formulae 
for the perturbative expansion of the tensor C"^^j given 
by Eq. (|3.3p . The tensor 6"^^^^ is also expanded in the 
same form as Eq. (|2.1ip . The first-order perturbations of 
C^^j have the well-known form(23] 

(I'C^,, = ^iaK^ ~ l^'Kb =: H^, ^ [h] , (B4) 

where if^^ [A\ is defined by Eq. p. lip for an arbitrary 
tensor field Aab defined on the background spacetime 
A^o- In terms of the tensor field H^^^'^ defined by p. lip 
the second-order perturbation ^^^C'^ab the tensor field 
C'afc is given by 

^■'^C^,, = H,,m]~2h^''HaM[h]. (B5) 

The Riemann curvature (|3.4p on the physical spacetime 
A^A is also expanded in the form ()2.1ip : 

n d n d , d , ,lx2(2)p d 

^abc —■ ^abc + ^ ^abc ' ' 2 '^'"^ 

+0(A3). (B6) 

The first- and the second-order perturbation of the Rie- 
mann curvature are given by 

= -2V[„(i)C'^,]„ (B7) 
= ~2ViJ'^C% + A^'^C\[J'^C% (B8) 

Substituting Eqs. (HH and into Eqs. (HZ]) and (|B8]) . 
we obtain the perturbative form of the Riemann curva- 
ture in terms of the variables defined by Eq. (|3.1ip and 
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= -2V[aH^/[h], (B9) 

+Ah''''y[aHb]ce [h] . (BIO) 

To write down the pcrturbative curvatures (|B9p and 
(jBlOp in terms of the gauge invariant and variant vari- 
ables defined by Eqs. p.23p and p.32p . we first derive 
an expression for the tensor field Habc[h] in terms of the 
gauge invariant variables, and then, we derive a pertur- 
bative expression for the Riemann curvature. 

First, we consider the linear-order perturbation (IB9p of 
the Riemann curvature. Using the decomposition (I2.23P 
and the identity R^g^^c] = 0, we can easily derive the 
relation 

Habc [h] = Habc [H] + Va^bX, + i?^,, "X^, (Bll) 

where the variable Habc [H] is defined by Eqs. (|3.1f p and 
([XT^ with Aab = Tiab- Clearly, the variable i/^^ " [H] 
is gauge invariant. Taking the derivative and using the 
Bianchi identity y[aRbc]de = 0, we obtain 

Wi?,,/ = -2V[„H,^/[n] + £xRabJ- (B12) 

Similar but some cumbersome calculations yield 

(2)i?,,/=-2V[„i7,j/ [£] + 4i/[/'= [H] Hb^cc [H] 
+m/V[aH,^/ [H] 

+2£x^'Kb/ + {£y - £l) Rabc'- (B13) 

Equations (jB12p and (|Bf 3p have the same for as the de- 
composition formulae (|2.36p and (|2.37p . respectively. 

Contracting the indices h and d in Eqs. (jBf 2p and 
(|B13p of the perturbative Riemann curvature, we can di- 
rectly derive the formulae for the perturbative expansion 
of the Ricci curvature: expanding the Ricci curvature 

Rab =■■ Rab + A(l)i?ah + +^A2(2)i?^, + 0{\^), (B14) 

we obtain the first-order Ricci curvature as 

<-''^Rab = -2VyaH,^b"M + £xRab- (B15) 



and we also obtain the second-order Ricci curvature as 

^ -2V^aH,^b"['^\+^H^a'^[nHc]bdm 

+ 2£x^'^Rab+{£Y-£x)Rab- (B16) 

The scalar curvature on the physical spacetime M is 
given hy R — g°'^Rab- To obtain the perturbative form of 
the scalar curvature, we expand the R in the form p. lip , 
i.e., 

i?=:i? + A(i)i?+^A2(2)i? + 0(A3) (B17) 

and g"'''Rab is expanded through the Leibniz rule. Then, 
the perturbative formula for the scalar curvature at each 
order is derived from perturbative form of the inverse 
metric (|B3P and the Ricci curvature (jBlSP and (|B16p . 
Straightforward calculations lead to the expansion of the 
scalar curvature as 

(i^i? - -2V[„F,] [H] - Rab'H'''' + £xR, (B18) 

= -2V^aHb-^'^'' [C]+ R''\2HcaHb' - Cab) 

+AH^a'^ [H] \ m + 4H,^V[,i/,] - m 

+47^''V[,i7,j,'^[H] 

+2£x^^^R+ {£y - £%:)R- (B19) 

We also note that the expansion formulae (|B18p and 
(|B19p have the same for as the decomposition formulae 
((OS)) and (|07)) . respectively. 

Next, we consider the perturbative form of the Einstein 
tensor Gab ■= Rab — \gabR and we expand Gab as in the 
form (Eini): 

Gab Gab + \^^KGab) + \\^^''HGab) + 0{X').{B2Q) 

As in the case of the scalar curvature, straightforward 
calculations lead 



^^\Gab) 
^'\Gab) 



-2V[,i7rfj/ [H] + gaby[cHa{^'- m - -RUab + ^gabRcd'H''' + £xGab, 
-2V[,i7,], ^ [£] + 4i/[/^ [H] H,]bd [H] + m/V[aH^^, ^ [U] 
^ ^ -2V[,iI,]^'^ [C] + 2Rden/H'''^ - RdeC'''' + 4i/[/^ [H] H^^\ [U] 



(B21) 



:9ab 



m/V[,H^^ {U\ + 4H^'=V[,i/<i], {n\) + 2Uab^[cild\ + nabU^'Rcd - Mab 



+ 2£x^^\Gab)+{£Y- £\)Gab- 



(B22) 
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We note again that Eqs. (jB2ip and (|B22[) have the same 
form as the decomposition formulae (|2.36p and (|2.37p . 
respectively. 

The perturbative formulae for the perturbation of the 
Einstein tensor 

G^" = f^Gac (B23) 
is derived by the similar manner to the case of the per- 



turbations of the scalar curvature. Through these formu- 
lae summarized above, straightforward calculations leads 
Eqs. (|3^ - (|3T0|) . We have to note that to derive the 
formulae p. 91) with Eq. p.lOp . we have to consider the 
general relativistic gauge-invariant perturbation theory 
with two infinitesimal parameters which is developed in 
Refs. [sl, I9II, as commented in the main text. 
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